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FEDERAL  AVIATION  ADMINISTRATION 
SYSTEMS  RESEARCH  AND  DEVELOPMENT  SERVICE 
SPECTRUM  MANAGEMENT  STAFF 


STATEMENT  OF  MISSION 

The  mission  of  the  Spectrum  Management  Staff  is  to  assist  the  Department  of  State, 
Office  of  Telecommunications  Policy,  and  the  Federal  Communications  Commission  in 
assuring  the  FAA's  and  the  nation's  aviation  interests  with  sufficient  protected 
electromagnetic  telecommunications  resources  throughout  the  world  to  provide  for  the  safe 
conduct  of  aeronautical  flight  by  fostering  effective  and  efficient  use  of  a natural 
resource-the  electromagnetic  radio-frequency  spectrum. 

This  objective  is  achieved  through  the  following  services: 

• Planning  and  defending  the  acquisition  and  retention  of  sufficient  radio-frequency 
spectrum  to  support  the  aeronautical  interests  of  the  nation,  at  home  and  abroad,  and 
spectrum  standardization  for  the  world's  aviation  community. 

• Providing  research,  analysis,  engineering,  and  evaluation  in  the  development  of 
spectrum  related  policy,  planning,  standards,  criteria,  measurement  equipment,  and 
measurement  techniques. 

• Conducting  electromagnetic  compatibility  analyses  to  determine  intra/inter-system 
viability  and  design  parameters,  to  assure  certification  of  adequate  spectrum  to  support 
system  operational  use  and  projected  growth  patterns,  to  defend  the  aeronautical 
services  spectrum  from  encroachment  by  others,  and  to  provide  for  the  efficient  use  of 
the  aeronautical  spectrum. 

•Developing  automated  frequency-selection  computer  programs/routines  to  provide 
frequency  planning,  frequency  assignment,  and  spectrum  analysis  capabilities  in  the 
spectrum  supporting  the  National  Airspace  System. 

• Providing  spectrum  management  consultation,  assistance,  and  guidance  to  all  aviation 
interests,  users,  and  providers  of  equipment  and  services,  both  national  and 
international. 
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PREFACE 
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EXECUTIVE  SUMMARY 


The  Electromagnetic  Compatibility  Analysis  Center  (ECAC)  has 
developed  an  analysis  model  (AVPAK)  to  determine  the  mutual  effects 
of  introducing  new  avionics  equipment  to  an  existing  airframe 
containing  operational  equipment.  This  was  accomplished  in  response 
to  a request  by  the  Federal  Aviation  Administration  (FAA) . The 
present  model  (AVPAK  3)  is  the  result  of  the  third  phase  of  work  in 

* this  area. 

The  analysis  of  the  mutual  effects  of  the  operation  of  equip- 

* ment  on  an  airframe  is  accomplished  by  predicting  the  expected  level 
of  interference  relative  to  the  degradation  threshold  of  each  receiver. 
Antennas  are  assumed  to  be  isotropic  and  may  be  collocated  on  the 
aircraft  or  on  a neighboring  aircraft.  Equipment  examined  for 
possible  interference  include  those  with  overlapping  or  immediately 
adjacent  operating  frequencies  and  also  those  with  harmonically 
related  operating  frequencies  for  which  inadequate  transmitter 
harmonic  attenuation  exists.  Nonlinear  effects  are  not  included 

in  the  analysis  and  must  be  dealt  with  manually.  To  insure  that 
only  far  field  coupling  conditions  are  considered,  only  those 
equipments  operating  above  30  MHz  are  treated. 

The  model  offers  the  option  of  either  a purely  deterministic 
calculation  or  a probabilistic  calculation  that  estimates  the  prob- 
ability of  interference.  The  body  of  the  airframe  is  modeled  as  a 
cylinder  of  finite  length  to  which  appropriate  conical  sections  may 
be  added.  The  calculation  takes  into  account  the  factors  of  air- 
frame curvature,  airfoil  obstructions,  and  bulkhead  obstruction. 


■ 


Two  data  bases  are  associated  with  the  model.  One  contains 
nominal  characteristics  of  commonly  used  avionics  equipment  which 
may  be  directly  called  upon  to  run  the  model.  The  model  contains 
a mode  which  allows  user-specified  data  to  be  used  in  running  the 
program  for  new  equipment  types  or  equipment  not  described  in  the 
data  base.  The  other  data  base  contains  median  values  and  standard 
deviations  of  selected  characteristics  for  equipments  grouped  accord- 
ing to  function.  One  set  of  characteristics,  based  on  median  values, 
is  then  used  to  represent  the  group.  This  information  is  used  in 
probabilistic  processing. 

As  part  of  this  task  the  existing  AVPAK  data  bases  were  expanded 
with  new  avionics  equipment  obtained  since  the  last  updating  in  1973. 

In  addition  to  interference  analysis,  the  model  calculates 
the  power  density  at  user-specified  points,  resulting  from  the 
operation  of  transmitters  located  on  an  aircraft.  These  points 
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EXECUTIVE  SUMMARY  (Continued) 


»ay  be  located  anywhere  on  the  airframe,  including  wing  pods,  or 
they  may  be  "raised"  from  the  airframe  (i.e.,  not  lying  on  the 
fuselage  skin),  including  locations  on  neighboring  aircraft.  The 
model  also  has  the  capability  of  calculating  a cumulative  power 
density  due  to  the  effects  of  more  than  one  transmitter. 
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The  portion  of  the  AVPAK  data  base  containing  equip- 
ment data  to  support  deterministic  analysis  of  EMC 
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of  Transportation. 

Frequency  Analysis  System  --  an  ECAC  analysis  program 
that  includes  the  capability  to  synthesize  transmitter 
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SECTION  1 
INTRODUCTION 


BACKGROUND 

In  1969,  the  Federal  Aviation  Administration  tasked  ECAC  to 
develop  an  automated  interference  analysis  capability  for  use  with 
avionics  equipment  in  aircraft.  The  capability  was  for  general 
application,  but  particularly  useful  when  a new  equipment  was  to  be 
introduced  to  an  existing  airframe  containing  operational  equip- 
ment. That  effort  was  completed  in  December  1970. 1 

In  1972,  FAA  requested  that  the  capability  be  extended,  notably 
to  allow  probabilistic  analyses.  Certain  other  improvements  were 
also  added,  and  additional  measured  data  was  obtained  for  further 
validation  of  the  antenna  coupling  portion  of  the  model.  The 
report  of  that  effort  was  published  in  1973. 2 

In  1974,  FAA  again  requested  that  ECAC  undertake  a task  for 
improving  the  model  and  incorporating  features  that  had  become 
available  since  publication  of  the  second  report.  A comparison 
of  this  report  with  References  1 and  2 will  reveal  much  repetitious 
material.  It  was  necessary  to  include  such  material  here  in  order 
to  satisfy  the  FAA  requirement  for  a single  document  describing 
the  capability.  Future  reference  to  the  previous  reports  for  an 
understanding  of  the  application  and  use  of  the  prediction  capability 
and  data  base  will  be  unnecessary. 

Throughout  this  report,  where  it  is  desired  to  distinguish 
between  the  early  model,  the  intermediate  model  and  this  latest 
model,  they  will  be  characterized  as  AVPAK  1,  AVPAK  2 and  AVPAK 
3,  respectively.  It  should  be  understood  that  AVPAK  1 and  AVPAK  2 
no  longer  exist.  AVPAK  3,  of  course,  has  all  the  capabilities  of 
the  first  two  versions  as  well  as  the  improvements  described  herein. 

OBJECTIVES 

The  objectives  of  this  project  were  to: 

1.  Incorporate  a power  density  model  in  AVPAK. 

2.  Incorporate  an  airfoil -shading  model. 


1 Morgan,  G.,  Avionic*  Interference  Prediction  Model,  ESD-TR-70-286, 
December  1970.  (FAA  Report  No.  FAA-RD-71-10.) 

2Friske,  L. , An  Extended  Avionics  Interference  Prediction  Model, 
ECAC- PR-73 -002,  June  1973.  (FAA  Report  No.  FAA-RD-73-9.) 
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3.  Develop  a model  for  predicting  coupling  between  two 
antennas  whose  locations  are  raised  from  the  airframe. 

4.  Expand  the  avionics-equipment  data  base. 

5.  Refine  and  regroup  equipment  data  to  arrive  at  smaller 
variances  for  use  in  probabilistic  analyses. 

6.  Validate  the  above  additions  to  the  model. 

7 . Document  the  prediction  model  and  data  base  to  provide 
a single  technical  report  describing  the  capability. 

APPROACH 

The  following  steps  were  taken  to  accomplish  the  specific 
objectives  of  the  project. 

1.  Certain  techniques  for  computing  air-foil  shading 
were  known  to  exist.  These  were  analyzed  and  a validation  was 
undertaken.  They  were  found  to  be  unsatisfactory  and  an  empirical 
model  was  developed. 

2.  A technique  for  computing  power  density  was  incorporated 
into  AVPAK . 

3.  Additional  manufacturers'  technical  manuals  were  obtained. 
Data  were  extracted  from  the  manuals  for  inclusion  in  the  AVPAK  data 
bases.  New  classifications  were  established,  leading  to  smaller 
variances  in  the  probabilistic  data  base. 

4.  Additional  coupling  loss  data  were  obtained  and  addi- 
tional verification  of  the  coupling  model  was  performed. 
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SECTION  2 
MODEL  EXTENSION 


GENERAL 

The  purpose  of  this  section  is  to  discuss  changes  made  to  the 
analysis  and  data  base  capabilities,  to  report  the  options  that 
were  open,  and  to  justify  the  selections  and  decisions  finally 
made.  Section  3 of  this  report  describes  the  entire  capability 
as  it  now  exists. 

DEVELOPMENT  OF  AN  AIRFOIL  OBSTRUCTION  LOSS 

The  initial  plan  for  satisfying  the  airfoil  obstruction- loss 
task  was  to  investigate  existing  obstruction  models  and  theories 
and  determine  if  their  use  in  the  AVPAK  model  was  feasible. 

The  first  area  investigated  involved  work  done  by  the  McDonnell 
Douglas  Corporation  and  the  Grumman  Aerospace  Corporation  in  two 
separate  ventures  where  wing  obstruction  attenuation  was  determined 
through  diffraction  effects.  Because  of  the  necessity  to  limit 
the  detail  associated  with  input  parameters,  it  was  not  feasible 
to  incorporate  either  of  these  models  in  the  AVPAK  model. 

An  attempt  was  made  to  model  an  airfoil  as  a knife-edge  obstruc- 
tion. When  this  yielded  inaccurate  results,  a decision  was  made 
to  develop  a model  to  account  for  obstructed  paths  — one  which 
could  be  easily  assimilated  into  the  AVPAK  model.  Each  of  these 
developments  in  the  obstruction-loss  task  is  discussed  below. 

McDonnell  Douglas  Model  (ATACAP) 

An  antenna -to -antenna  compatibility  analysis  program  (ATACAP) 
was  developed  by  the  McDonnell  Douglas  Corporation  under  USAF 
contract  as  part  of  a larger  intra-vehicle  electromagnetic  compat- 
ibility analysis  program. 

ATACAP,  like  the  Avionics  Interference  Prediction  Model, contains 
not  only  antenna-coupling  routines  but  receiver  and  transmitter 
synthesis  models.  Only  the  antenna  coupling  was  of  interest  here 
since  ATACAP  includes  provision  for  predicting  the  additional  loss 
introduced  by  wings.  This  additional  loss  is  based  on  the  Geometrical 


3Bogdanor,  J.  L.,  Siegal,  M.  D.,  Weinstock,  G.  L.,  Intra-Vehicle 
Electromagnetic  Compatibility  Analysis,  Part  I,  McDonnell  Air- 
craft Company,  McDonnell  Douglas  Corporation,  TR  AFAL-TR-71  155 
PTI,  January  1972. 
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Theory  of  Diffraction  (which  is  explained  in  some  detail  below  and 
in  APPENDIX  A) . 

Documentation  and  a program  listing  for  ATACAP  were  obtained 
and  evaluated  for  applicability  to  the  requirements  of  FAA.  The 
evaluation  revealed  apparent  irregularities  that  could  not  be 
accounted  for  either  by  analysis  or  consultation  with  McDonnell- 
Douglas . For  that  reason,  the  ATACAP  effort  was  terminated  and 
an  attempt  was  made  to  use  basic  work  on  the  Geometric  Theory  of 
Diffraction. 

Geometric  Theory  of  Diffraction 

The  Geometric  Theory  of  Diffraction  is  an  extension  of  geo- 
metrical optical  theory  which  has  been  widely  used  to  solve  antenna 
problems;  e.g.,  the  determination  of  antenna  patterns  for  multi- 
element arrays,  high  gain  reflectors  and  dielectric  lenses. 

Geometric  optics  theory,  however,  fails  to  account  for  dif- 
fraction which  occurs  when  the  incident  ray  hits  edges,  corners, 
apexes  of  surfaces,  or  is  tangent  to  a smooth  surface.  GTD  was 
developed  to  account  for  these  optical  phenomena  of  diffraction.4 

An  exposition  of  the  ECAC  inquiry  into  GTD  for  modeling 
obstuctions  between  aircraft  antennas  appears  in  APPENDIX  A. 

It  was  determined  that  the  geometric  detail  necessary  to  use  such 
a model  would  be  available  only  after  precision  measurements  of 
aircraft  geometry.  Because  of  this  limitation  GTD  was  rejected 
as  an  acceptable  obstruction-modeling  technique. 

Knife  Edge  Diffraction 

The  classic  Bullington  equation,  used  by  the  model  to  cal- 
culate transmission  loss  over  a bulkhead  obstruction,  was  con- 
sidered initially  for  calculating  the  loss  around  an  airfoil  obstruc 
tion.  However,  it  was  found  that  this  technique  does  not  yield 
verifiable  results.  Therefore,  it  was  abandoned  for  calculating 
airfoil  obstruction  losses. 

Empirical  Airfoil  Obstruction  Model 

Since  none  of  the  investigated  methods  yielded  the  desired 
results  in  terms  of  accuracy  in  predicting  obstructed-path  coupling 


4Sacks,  L.H.,  The  Geometrical  Theory  of  Diffraction  Applied  to 
Aircraft  Antenna  Isolation  Determination,  Parts  I and  II,  No. 
RF-67-10,  Grumman  Aerospace  Corp. , May  1967. 
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losses,  an  empirical  method  to  compute  airfoil  obstruction  loss 
was  developed  at  ECAC.  Through  inspection  of  measured  path 
losses  where  airfoil  obstructions  existed,  it  was  determined  that 
obstruction  loss  could  be  approximated  by  calculating  the  free- 
space  loss  around  the  airfoil  and  adding  a curvature  factor  (see 
Figure  1).  The  free-space  loss  is  obtained  from  the  following 
equation: 

Lps  = 20  log  D + 20  log  f -38  (1) 

where 

LpS  * free-space  loss  in  dB 
D = minimum  distance  around  airfoil,  in  feet 
f * transmitter  frequency,  in  MHz. 

The  curvature  loss  is  obtained  from  the  curvature  loss  function 
F(y),  which  is  discussed  in  Section  3.  For  an  airfoil  obstruction, 
in  calculating  F(y),  the  length  of  the  ray  path  between  the  two 
antennas  is  the  minimum-distance  path  between  the  antennas,  (through 
the  airfoil). 

The  chords  of  all  airfoils  are  assumed  to  be  parallel  to  the 
Z-axis.  If  the  airfoil-fuselage  chord  does  not  lie  parallel  to 
the  Z-axis,  the  chord  will  be  set  parallel  to  the  Z-axis  in  the 
program,  with  reference  to  the  front  airfoil-fuselage  intersection 
point.  An  airfoil  obstruction  loss  will  be  calculated  if  antennas 
are  located  in  the  following  areas  (see  Figure  2).  Both  antennas 
must  be  located  within  ±20  degrees  of  the  top  or  bottom  fuselage 
centerline  and  lie  on  the  fuselage.  They  also  must  be  located  so 
that  there  is  full  airfoil  obstruction.  This  is  determined  by 
comparing  the  Z-axis  coordinates  of  the  two  antennas  with  those 
of  the  forward  and  aft  wing-fuselage  intersection  points.  Let 
(Xj,  y^,  Zj)  and  (x2,  y2  > z2)  t*'e  coordinates  of  the  two  antennas 

and  (x3,  yy  z^)  and  (x4,  y4,  z4)  the  coordinates  of  the  forward 

airfoil-fuselage  intersection  point  and  aft  airfoil-fuselage  inter- 
section point,  respectively.  Airfoil  shading  is  calculated  if: 

1)  (z3  - e)  < Zj  < (z4  + e)  and  z3  < z2  < z4  (2a) 

or 

(z3  - e)  < *2  < (z4  ♦ c)  and  z3  < Zj  < z4  (2b) 


2) 
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where 


e - (-1)  (z4  - z3). 

This  technique  provides  accurate  wing-obstruction  calculations 
using  a minimum  number  of  readily  available  user  inputs. 

RECLASSIFICATION  TO  REDUCE  STANDARD  DEVIATIONS 

One  of  the  objectives  of  the  previous  ECAC  work  on  the  Avionics 
Interference  Prediction  Model  entailed  the  creation  of  a functional 
data  base  comprised  of  representative  equipments.  The  parameters 
of  these  representative  equipments  were  expressed  in  terms  of  mean 
values  and  their  associated  standard  deviations.  A functional 
representation  is  used  in  the  probabilistic  mode  of  analysis. 

Since  large  standard  deviations  were  associated  with  some  of  the 
functional  (representative)  equipments,  the  task  was  undertaken 
to  reclassify  the  data  base  elements  such  that  probabilistic 
analyses  could  be  performed  with  a higher  degree  of  confidence 
in  the  predicted  interference  levels  (i.e.,  smaller  standard  devia- 
tions) . 

It  was  suggested  by  the  FAA  that  either  the  FAA  Technical 
Standard  Order  (TSO)  or  ARINC  Equipment  Characteristics  (ARINC 
EC)  be  considered  as  a means  of  classifying  the  equipment.  TSO's 
contain  minimum  performance  and  quality-control  standards  for 
specified  materials,  parts  or  appliances  used  on  civilian  aircraft. 
They  are  issued  by  the  Federal  Aviation  Administration  (FAA) , and 
the  performance  standards  in  each  TSO  are  those  that  the  FAA  finds 
necessary  to  ensure  that  the  equipment  concerned  will  operate 
satisfactorily  or  will  accomplish  satisfactorily  its  intended 
purpose  under  specified  conditions.5 

An  ARINC  EC  is  used  to  indicate  to  prospective  manufacturers 
of  airline  electronic  equipment  the  considered  opinion  of  the  air- 
line industry  concerning  requisites  for  new  equipment.  They  also 
serve  the  purpose  of  standardizing,  to  the  maximum  possible  extent, 
those  physical  and  electrical  characteristics  which  affect  inter- 
changeability of  equipment.6  Both  the  TSO  and  ARINC  EC  were 
examined  to  see  whether  samples  of  equipments  certified  to  meet 


5Federal  Aviation  Regulations,  Part  37,  Technical  Standard  Order 
Authorization,  May  1974,  Department  of  Transportation,  Federal 
Aviation  Administration. 

6ARINC  Characteristic  #579-1.  Aeronautical  Radio,  Inc.,  2551 
Riva  Rd,  Annapolis,  MD  21401,  5 February  1971. 
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the  requirements  of  either  (or  both)  of  them  were  characterized 
by  smaller  (or  larger)  standard  deviations.  No  such  correlation 
could  be  found.  In  some  cases,  for  example,  sharper  selectivity 
characteristics  are  associated  with  communications  receivers  that 
are  not  certified  to  meet  either  a TSO  or  ARINC  EC  than  with  equip- 
ment so  certified. 

Nominal  characteristics  provided  by  manufacturers  do  not 
always  indicate  the  true  characteristics  of  equipment.  That  is, 
nominal  values  are  cited  for  selectivity,  for  example,  that  meet 
some  specification  requirement  when  in  fact  the  selectivity  may 
be  far  superior  to  that  representation.  This  practice  would, 
of  course,  make  variances  quite  wide  when  some  values  are  based 
upon  these  nominal  characteristics  and  others  are  actual  values. 

Other  approaches  were  taken  to  the  reclassification  problem. 
For  example,  equipments  were  ranked  by  standard  deviation  of 
pertinent  characteristics,  and  then  they  were  examined  for  other 
common  or  similar  characteristics.  It  was  decided  that  there  was 
no  nominally  stated  feature  upon  which  to  base  classification. 

Therefore,  equipment  types  associated  with  large  standard 
deviations  (viz.  VHF  communications  receivers  and  VOR/Localizer 
receivers)  were  subclassified  according  to  the  highest  recorded 
spurious-response  levels  in  each  category.  Three  subclasses 
were  used,  as  follows: 

Subclass  Relative  Spurious  Response  Level 

1 40-89  dB  down 

2 90  - 109  dB  down 

3 110  - 140  dB  down 

POWER  DENSITY  MODEL 


An  earlier  ECAC  effort  was  involved  with  the  evaluation  of 
the  electromagnetic  environment  at  different  points  on  an  air- 
craft. This  effort  developed  the  capability  of  calculating 
power  densities  resulting  from  transmitters  located  on  the  air- 
craft as  well  as  from  transmitters  located  on  neighboring  air- 
craft. This  model  has  been  incorporated  into  AVPAK  3.  Details 
on  the  power  density  model  are  contained  in  the  model  description 
section  of  this  report  (see  pg.  43). 
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SECTION  3 
MODEL  DESCRIPTION 


GENERAL 


In  determining  possible  interference  cases,  the  AVPAK  3 model 
utilizes  two  major  analyses.  One  calculates  the  path  loss  between 
isotropic  radiators  located  on  an  airframe,  raised  from  an  air- 
frame, or  on  neighboring  aircraft.  An  antenna  gain  subroutine 
provides  gain  in  the  direction  of  coupling  for  certain  types  of 
aperture  antennas  while  non-aperture  antennas  are  assigned  a nominal 
gain  value.  The  other  major  analysis  determines  the  rejection 
offered  by  a receiver  to  an  undesired  signal.  This  is  accomplished 
by  integrating  over  the  areas  of  frequency  overlap  between  the 
spectral  emission  of  a chosen  transmitter  and  the  receiver  response 
characteristic  of  a chosen  receiver.  For  equipment  such  as  ATC 
transceivers  that  operate  over  a range  of  frequencies,  the  entire 
operating  range  is  examined;  in  this  manner  "worst  case"  values 
for  rejection  are  determined. 

Also  included  in  the  AVPAK  model  is  the  capability  to  calculate 
power  density  at  receivers  or  user-specified  points  due  to  a number 
of  transmitters  (up  to  50) . 

Nonlinear  effects  such  as  intermodulation,  cross-modulation, 
receiver  desensitization  and  specific  spurious  responses  are  not 
included  in  AVPAK  3 and  must  be  treated  manually. 

In  creating  the  AVPAK  3 model,  the  desired  new  features  and 
capabilities  were  added  to  the  AVPAK  2 format.  The  specific 
capabilities  of  the  earlier  AVPAK  models  and  the  new  features 
included  in  the  AVPAK  3 model  are  discussed  below. 

BASIC  INTERFERENCE  EXPRESSION 

The  analysis  of  the  mutual  effects  of  the  operation  of  equip- 
ment on  an  airframe  is  accomplished  by  predicting  the  expected 
level  of  interference  relative  to  the  degradation  threshold  of 
each  receiver. 

The  undesired  interfering  power  at  the  input  terminals  of  a 
potential  victim  receiver  can  be  calculated  with  the  logarithmic 
form  of  the  one-way  system  loss  equation,  which  is: 

pi  “ pt  * «r  * gr  - lp  - Ls 
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where 

Pj  * the  input  interfering  power  in  dBm 

P.j.  = the  transmitter  output  power  in  dBm 

Gy  ■ the  effective  gain  of  the  transmitting  antenna  in 
the  direction  of  the  receiving  antenna  in  dBi 

Gr  * the  effective  gain  of  the  receiving  antenna  in  the 
direction  of  the  transmitting  antenna  in  dBi 

Lp  3 the  basic  transmission  loss  between  isotropic  radi- 
ators in  dB 

Lg  = the  combined  system  losses  associated  with  the  trans- 
mitter and  receiver  due  to  transmission  links,  coupling 
devices,  and  external  RF  filtering, in  dB. 

Except  for  coupling  mismatches,  the  system  losses  can  usually 
be  neglected.  It  is  assumed  herein  that  such  losses  are  negligible. 
Therefore,  the  Lg  term  is  dropped  from  further  consideration  at 

this  time,  but  the  capability  to  include  such  a factor  has  been 
retained  for  future  application. 

The  minimum  level  of  an  input  desired  signal  required  to  pro- 
duce a standard  response  is  defined  as  the  receiver  sensitivity, 

Rg.  The  level  of  an  interfering  signal  required  to  produce  the 

same  standard  response  in  the  receiver  may  be  different.  The  two 
factors  which  can  cause  this  difference  are  the  different  modulation 
characteristics  of  the  interfering  and  desired  signals,  and  the 
rejection  to  an  off-tune  interfering  signal  offered  by  the  selectivity 
characteristic  of  the  receiver.  These  factors  can  be  considered 
together  as  the  total  receiver  rejection  factor  (L^) . 

There  is  an  equivalent  on-tune  input  signal  poweT  that  would 
produce  the  same  response  at  the  output  as  the  interfering  signal. 

It  is  related  to  the  interfering  signal  by: 


where 

Pie  * the  equivalent  on-tune  input  signal  power  that 

produces  the  same  output  response  as  the  interfering 
signal.  (If  P^e  * Rs,  the  output  would  be  a stand- 
ard response.) 
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By  substitution  in  Equation  3: 

Pie  = PT  + GT  + GR  " LP  ' Lf 

The  level  of  degradation  caused  by  an  interfering  signal  can 
be  evaluated  by  comparing  the  resulting  signal-to-interference 
ratio  to  the  required  threshold  S/I  ratio.  If  it  is  assumed  that 
the  input  desired  signal  is  at  the  receiver  sensitivity  level, 
the  relative  degradation  level  is: 

Pld  *»'  Pie  - Rs  * <S/I>T  (6) 

where 

Pld  = t*ie  degradation  level  relative  to  the  threshold  of 
the  degradation  in  dB 

Rs  = the  sensitivity  of  the  receiver  in  dBm 

(S/I)T  = the  threshold  input  signal-to-interference  ratio, 
in  dB,  required  to  prevent  degradation. 

If  the  desired  signal  is  at  a level  other  than  R , the  actual  level 

s 

may  be  substituted  for  Rs>  If  Equations  5 and  6 are  combined: 

Pld  = PT  + GT  + GR  ‘ LP  ‘ Lf  + (S/I>T  - Rs  (7) 

where  all  of  the  terms  have  been  defined. 

If  Pj^  is  greater  than  zero,  degradation  is  expected  to  occur; 

conversely,  if  P ^ is  less  than  zero,  degradation  is  not  expected 

and  no  further  consideration  need  be  given  to  interactions  between 
the  particular  transmitter  and  receiver  involved.  Each  of  the 
terms  in  Equation  7 is  discussed  below. 

Transmitted  Power,  PT  (dBm) 


This  parameter  is  a required  input  to  the  program  and  represents 
the  average  output  power  for  communications  transmitters  and  the 
peak  output  power  for  pulsed  transmitters.  This  information  can 
be  obtained  from  either  of  the  AVPAK  3 data  bases  (AVBASE , or  AVFILE) 
or  it  may  be  provided  by  the  user. 
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Antenna  Gains,  G^,,  GR  (dBi) 


AVPAK  3 has  the  capability  of  utilizing  either  user-supplied 
antenna  gains  or  calculated  gains  determined  by  automated  antenna 
gain  subroutines  which,  for  different  arrays,  determine  the  ex- 
pected gain  to  be  realized  along  the  propagation  path  between  the 
antennas  on  the  airframe.  For  aperture  antennas,  gain  routines 
were  developed  for  three  basic  types:  horn,  circular  aperture 

and  rectangular  array.  These  routines  calculate  far-field  gain 
using  aperture  dimensions,  transmitter  frequency,  and  coupling 
angles  referenced  to  the  main  beam  axis.  By  calculating  receiver 
antenna  gain  at  the  transmittex  frequency  (either  fundamental  or 
harmonic),  out-of-band  antenna  responses  are  treated  for  aperture 
antennas.  For  non-aperture  antennas,  a nominal  (2  dBi)  gain  is 
assigned. 

The  following  material  briefly  describes  some  of  the  considera- 
tions in  deriving  antenna  gain  values,  and  the  equations  used. 

The  main  beam  gain  of  a circular  aperture  antenna  is  calculated 


in  the  program 

as: 

GMB  = 

20  log  D + 20  log  F - 52.6 

(8) 

where 

r 

gmb  ' 

main  beam  gain  in  dBi 

0 . 

diameter  of  aperture  in  feet 

1 

frequency  in  MHz. 

Sidelobe  gain  is  determined  by  using  a Bessel  envelope  calculation. 

The  main  beam  gain  for  a rectangular  aperture  antenna  is  calcu- 
lated from: 

Gmb  = 10  log  (LX  x LY)  + 20  log  F - 51.6  (9) 

where 

LX  and  LY  = aperture  dimensions  in  feet 
F = frequency  in  MHz. 

Sidelobe  gain  is  determined  from  a sinX/X  envelope  calculation. 7 


7Kraus,  J.  D.,  Antennas,  McGraw-Hill,  New  York,  1950. 
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The  main  beam  gain  for  an  optimum  horn  is  calculated  from: 
G^  = 8.08  + 10  log  (LA  x LB)/A2  (10) 

where 


LA  and  LB 


aperture  dimensions  in  inches 


\ = wavelength  in  inches. 


Sidelobe  levels  for  an  optimum  horn  are  determined  from  look- 
up tables  in  the  program.  The  method  for  computing  these  levels 
was  obtained  from  Jasik.8 


For  all  antenna  types,  backlobe  gain  is  constant  and  equal 
to  the  value  of  the  gain  envelope  at  the  90°  pattern  point.  All 
gains  are  idealized  far-field  values.  Pattern  perturbations 
which  may  be  caused  by  the  airframe  are  not  considered. 

Polarization  mismatches  are  applied  only  to  cross-polarized 
non-aperture  antennas  (i.e.,  for  vertical-to-horizontal  polarization 
or  vice  versa) . Polarization  mismatches  are  not  applied  to  aperture 
antennas  in  the  program,  since  this  correction  is  appropriate  for 
main  beam  coupling  only,  a condition  which  rarely  occurs.  When 
appropriate,  up  to  20  dB  of  loss  can  be  manually  applied  to  the 
results  for  cross-polarized  main  beam-coupled  antennas.  Mismatches 
for  circularly  or  linearly  (45°)  polarized  antennas  to  vertically 
or  horizontally  polarized  antennas,  (or  vice  versa)  are  not  treated 
in  the  program.  When  necessary  a 3-dB  loss  may  be  manually  applied 
to  the  results  of  such  a configuration. 

Cosite  Path  Loss,  Lp 


The  path  loss  between  isotropic  radiators  on  an  airframe  is 
calculated  using  the  technique  reported  by  Hasserjian  and  Ishimaru9 
and  extended  by  Khan,  et  al.10  These  efforts  have  shown  that  the 
path  loss  along  a conducting  curved  surface  can  be  calculated  by: 


8Jasik,  H.,  ed.,  Antenna  Engineering  Handbook,  McGraw-Hill,  New 
York,  1961. 

9Hasserjian,  G.,  and  Ishimaru,  A.,  Excitation  of  a Conducting 
Cylindrical  Surface  of  Large  Radius  of  Curvature,  IRE  Transactions 
on  Antennas  and  Propagation,  Vol  AP-10,  May  1962. 

1 °Khan,  P.  J.,  et  al.  Derivation  of  Aerospace  Antenna  Coupling 
Factor  Interference  Prediction  Techniques,  Colley  Electronics 
Laboratory,  University  of  Michigan,  1964. 
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LPC  = LPF  F(y)  (ID 

where 

Lpp  = the  path  loss  along  the  curved  surface  of  the  air- 
frame 

Lpp  * the  path  loss  if  the  same  surface  were  flattened 
into  a plane 

F(y)  = the  loss  factor  due  to  the  curvature  of  the  surface 
i.e.,  the  curvature  factor. 

The  curvature  factor  F(y)  was  expressed  by  Hasserjian, 
(Reference  9)  in  the  form  of  two  different  infinite  series  as  an 
approximate  solution  to  Maxwell's  equation  with  boundary  condi- 
tions. One  series  expression  was  derived  for  "large  values"  of 
the  parameter  y and  the  other  series  was  derived  for  "small  values" 
of  y.  The  numerical  evaluation  of  these  series  for  the  magnitude 
of  F (y)  versus  y in  decibels  was  done  at  the  University  of  Michigan 
(Reference  10),  and  is  presented  as  a tabulatipn  of  points  of  the 
graph  of  F (y)  versus  y in  Figure  3.  This  tabulation  is  used  with 
an  interpolation  subroutine  and  the  formula  for  y. 

The  independent  variable,  y,  of  this  function  is  in  itself  a 
function  of  the  "ray  path"  of  length  R,  and  of  curvature  p..  The 
function  y is:  1 

y = (R1)3/2/p1  (12) 

where 

Rj  = the  length  of  the  curved  ray  path  as  normalized  to 
the  wave  number,  k (where  k = 2tr/X,  and  X is  the 
wavelength).  That  is,  Rj  * k Dj,  where  D]L  is  the 

length  of  the  curved  ray  path 

Pj  = the  curvature  of  the  ray  path  as  normalized  by  the 
wave  number. 


The  value  of  y for  a cylinder  is  determined  from  the  following 
formula: 
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w 

where 

a * the  radius  of  the  cylinder 

AZ  » the  distance  between  the  antennas  along  the  axis 
of  the  cylinder 

| ♦ * the  angle  in  radians  between  the  antennas  on  a plane, 

• defined  by  the  two  antennas  and  the  center  of  the  air- 

frame. 


The  path  loss  between  the  antennas  if  the  surface  is  assumed 
to  be  planar  is  calculated  using  the  free-space  formula: 


lpf 

20  log  f ♦ 20  log  Dx  - 38 

(14) 

where 

lpf 

the  free-space  loss  between 
dB 

isotropic  radiators  in 

f * 

the  transmitted  frequency  in  MHz 

D1  = 

the  ray  path  distance  along 

the  surface  between  the 

antennas  in  feet. 

The  distance,  D^,  for  a cylinder  is: 

D1  * [(AZ)2  + (ad.)2]*4  (14a) 

where  all  terms  have  been  previously  defined. 

One  of  the.  restrictions  in  the  technical  development  of  the 
curvature  factor  is  a requirement  that  the  curvature  along  the 
ray  path  between  the  antennas  remain  constant.  If  the  antennas 
lie  on  a conical  surface,  the  requirement  is  not  completely 
satisfied.  In  practical  airframes,  however,  it  can  be  shown  that 
the  ray  path  length  can  be  calculated  with  a high  degree  of  accuracy 
by  treating  the  cone  as  a modified  cylinder  with  a radius  equal 
to  the  geometric  mean  of  the  radii  of  the  cone  at  the  locations 
of  the  antennas.  The  limitation  of  this  approximation  is  that  the 
apex  angle  of  the  cone  can  not  exceed  20  degrees. 

It  can  be  shown  that  the  curvature  factor  between  two  antennas 
on  a cone  lies  between  the  two  F(y)'s  which  would  be  calculated  if 
the  cone  was  replaced  by  two  different  cylinders  having  radii  equal 
to  the  cone  radii  at  each  of  the  two  antenna  locations. 
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The  remaining  restrictions  which  affect  the  application  of  this 
technique  include  geometrical  limitations  which  insure  that  the 
respective  antennas  do  not  lie  within  each  other's  Fresnel  (near- 
field)  region.  These  geometrical  restrictions  place  a lower  limit 
on  the  frequencies  at  which  the  coupling  loss  can  be  calculated. 

For  example,  if  an  HF  wire  antenna,  the  resulting  radiating  high 
frequency  wavelengths,  and  an  airframe  length  arc  all  of  comparable 
magnitudes,  then  the  entire  airframe  can  be  expected  to  be  a part 
of  the  HF  antenna  system.  Any  consideration  of  the  coupling  loss 
to  be  expected  along  the  airframe,  consequently,  becomes  an  intra- 
antenna (near-field)  consideration.  Thus,  this  technique  cannot 
be  applied  to  HF  systems. 

In  addition  to  the  parameters  discussed  above,  other  factors 
must  be  considered  for  certain  paths  between  antennas  on  an  air- 
frame. When  one  antenna  is  located  in  front  of  the  metal  nose 
bulkhead  and  the  other  is  behind  the  bulkhead,  a knife-edge  dif- 
fraction loss  can  be  expected  due  to  the  obstruction  created  by  the 
bulkhead.  Bullington  presents  a nomograph  which  can  be  used  to 
calculate  these  losses  (Reference  8,  Chapter  33).  This  nomograph 
has  been  automated  in  equation  form  and  this  additional  loss 
factor  is  automatically  included  when  the  transmission  path  crosses 
the  nose  bulkhead.  The  equation  used  is: 


Lk  = 

10  log  (h2f/20d) 

(15) 

where 

Lk  * 

the  knife-edge  diffraction  loss,  in  dB 

h «= 

the  height  of  the  obstruction  above  the 
line  path,  in  feet 

straight- 

f * 

the  transmitted  frequency,  in  megahertz 

d » 

the  distance  between  the  bulkhead  and  the  antenna 
nearest  the  bulkhead,  in  feet. 

Another  situation  which  must  be  considered  occurs  when  the 
path  between  antennas  is  obstructed  by  an  airfoil.  A quasi-empirical 
equation  has  been  developed  at  ECAC  which  calculates  the  coupling 
loss  of  an  airfoil-obstructed  path  by  determining  the  minimum  free- 
space  loss  around  the  obstruction  and  adding  a curvature  loss  over 
the  minimum-distance  path  between  antennas  (through  the  airfoil). 

L - 20  log  f ♦ 20  log  D - 38  ♦ F(y)  (16) 

W 

where 

Lg  * wing  obstruction  loss  in  dB 
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f = transmitted  frequency  in  MHz 
D = minimum  distance  around  wing  in  feet 


F(y)  = previously  discussed  curvature  factor  whose  ray 
path  distance  neglects  the  presence  of  a wing 
(see  Section  2 for  additional  development  of 
wing-obstructed  coupling  loss). 


AVPAK  3 can  also  handle  transmitters  or  receivers  (or  points, 
if  a power  density  calculation  is  desired)  that  are  raised  from 
the  fuselage;  Taised  is  defined  as  not  lying  on  the  fuselage  skin. 
AVPAK  3 will  calculate  coupling  loss  for  the  following  configurations: 

1.  Both  antennas  lie  on  the  fuselage. 

2.  Both  antennas  are  raised  from  the  fuselage. 

3.  Antenna  1 is  raised  from  the  fuselage,  antenna  2 lies 
on  the  fuselage. 

4.  The  path  between  antenna  1 and  antenna  2 is  obstructed 
by  an  airfoil. 

5.  Antenna  1 lies  on  the  fuselage,  antenna  2 lies  on  a 

wing  pod. 

6.  Antenna  1 is  raised  from  the  fuselage,  antenna  2 lies 
on  a wing  pod. 


Intersite  Path  Loss 


For  intersite  analysis  the  path  loss  is  calculated  using  the 
smooth  curve-smooth  earth  (SCSE)  model  which  is  part  of  ECAC's 
master  propagation  calculation  system.  The  entire  airframe  is  assumed 
to  be  a radiator;  hence,  the  individual  bulkhead,  curvature,  and  wing- 
shading  losses  associated  with  the  cosite  path  are  assumed  to  be 
negligible.  The  SCSE  model11  assumes  a smooth  spherical  earth  with 
a 4/3  radius  of  earth  curvature.  The  model  is  valid  for  frequencies 
in  the  SO  MHz- to- 12  GHz  range  and  is  applicable  for  antenna  heights 
of  great  altitudes. 

SCSE  computes  the  basic  median  far- field  path  loss  between 
isotropic  radiators,  using  effective  antenna  heights  at  the  two 
site  locations.  The  model  automatically  selects  the  appropriate 
propagation  region  to  be  analyzed:  reflection,  intermediate,  or 

diffraction.  If  antenna  heights  are  greater  than  3,000  feet,  which 
is  usually  the  case  for  an  aircraft  in  flight,  a preprocessor  of 


^Leggett,  Robert  and  Madison,  James,  Propagation  User's  Manual, 
ECAC-UM-74- 001,  ECAC,  Annapolis,  MD,  July  1974. 
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the  model  calculates  a ray-trace  correction  along  with  the  effective 
earth  modeling.  Additional  information  about  the  SCSE  model  may 
be  found  in  Reference  11. 

In  an  intersite  analysis  it  is  possible  to  have  an  aircraft- 
to-satellite  situation,  where  the  aircraft  has  a flight  altitude 
of  less  than  100,000  feet  (approximately  19  statute  miles)  and 
the  satellite  an  altitude  greater  than  1,550,000  feet  (approximately 
294  statute  miles).  In  a case  such  as  this,  an  alternative  approach 
to  SCSE  is  used. 

The  geometry  of  the  aircraft-to-satellite  analysis  situation 
is  shown  in  Figure  4.  The  initial  points  and  distances  used  are: 

A.  location  of  aircraft  in  flight  (<  100,000  feet)  approx- 
imated to  be  on  the  earth's  surface  relative  to  the  altitude  of 

the  satellite. 

B.  location  of  the  satellite. 

C.  center  of  the  earth. 

E.  point  on  the  earth's  surface  on  the  line  BC. 

S.  ground  distance  (Great  Circle)  from  A to  E along  the 
surface  of  the  earth. 

CE.  approximate  radius  of  the  earth,'  (3963.0  statute 

miles). 

BE.  altitude  of  the  satellite  (>>  1,550,000  feet). 

For  the  aircraft-to-satellite  case,  BE  is  usually  much  greater 
than  1,550,000  feet.  Because  of  this,  the  altitude  of  the  aircraft 
in  flight  (while  in  fact  some  height  less  than  100,000  feet)  may 
be  approximated  as  0 feet.  If  the  altitude  of  the  aircraft  is 
assumed  to  be  0 feet,  then  the  aircraft  lies  on  the  curve  S along 
the  surface  of  the  earth. 

Let  BE  represent  the  altitude  of  the  satellite.  If  the  segment 
S is  less  than  or  equal  to  500  statute  miles  then  the  path  loss  is 
calculated  by: 

Lp  - 20  log  (f)  + 20  log  (BE)  +37.0  (17) 

where 

Lp  - path  loss  between  points  A and  B,  in  dB 

f » transmitter  frequency,  in  MHz 

BE  ■ altitude  of  the  satellite,  in  statute  miles, 
and  BE  - BA 
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Figure  4.  Aircraft -to-satellite  analysis  geometry.  A fsgnimt 

of  the  earth's  surface  between  A and  E is  represented 
by  the  curved  line  S.  Point  A is  the  aircraft  in 
flight,  which,  relative  to  the  satellite  at  point  B, 
can  be  approximated  to  lie  on  line  S. 
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This  is  the  free-space  loss  equation  for  a line-of-sight  case, 
which  is  valid  when  S is  less  than  or  equal  to  500  statute  miles. 

If  the  segment  S is  greater  than  500  statute  miles,  it  must 
be  determined  whether  the  two  sites  are  within  line-of-sight  of 
each  other.  This  particular  analysis  will  treat  only  line-of-sight 
cases.  A detailed  analysis  may  be  found  in  Reference  12.  If 
point  E is  found  to  be  beyond  line  of  sight  from  point  A,  the  run 
is  aborted.  If  not,  Lp  is  calculated  by  the  free-space  loss  equation 

with  a more  accurate  representation  of  the  distance  between  the 
sites  AE,  than  is  used  in  Equation  17  (Reference  11). 

Receiver  Rejection, 


The  rejection  offered  to  an  undesired  emission  by  a potential 
victim  receiver  is  calculated  by  a programmed  package  known  as 
the  Frequency  Analysis  System  (FAS).  The  following  explanation 
of  FAS  will  enable  the  reader  to  understand  the  operation  of  FAS 
and  will  prepare  him  for  the  more  detailed  description,  including 
the  mathematical  approach,  contained  in  the  documentation  of  the 
model. 13 

For  a given  transmitter-receiver  pair,  FAS  synthesizes  the 
receiver  response  characteristic  and  the  transmitter  spectral  emission 
characteristic  by  a series  on  line  segments  which  are  linear  on  a 
logarithmic  scale.  Each  of  these  synthesized  models  is  normalized 
to  unity  at  the  tuned  frequency  of  the  equipment  such  that  the 
characteristics  described  are  relative  to  the  performance  at  the 
tuned  frequency.  The  FAS  synthesis  does  not  include  nonlinear 
effects.  When  this  synthesis  is  complete,  FAS  determines  the  rejec- 
tion by  integrating  over  the  areas  of  frequency  overlap  between 
the  transmitted  emission  and  the  receiver  response.  Two  cases 
are  considered.  The  first  case  examined  involves  a calculation 
of  the  relative  energy  transfer  to  be  expected  due  to  the  emission 
sidebands  which  lie  within  the  passband  of  the  receiver.  The 
second  case  examined  is  the  expected  energy  transfer  resulting 
from  inadequate  receiver  selectivity  at  the  frequencies  within 
the  fundamental  emission  bandwidth  of  the  transmitter.  These 


12Haseltine,  R.,  Avionics  Interference  Prediction  Model  (AVPAK) , 
ECAC-TN-75-020,  ECAC,  Annapolis,  MD,  September  1975. 

13Cleaver,  R.  and  Bode,  T. , An  Algorithm  for  Calculating  Transmitter- 
Receiver  Frequency  Rejection  Loss,  ESD-TR-70-128,  ECAC,  Annapolis, 
MD,  1970. 
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two  cases  are  compared  and  the  worst  situation,  i.e.  the  least 
amount  of  rejection,  is  chosen  as  the  appropriate  rejection  for 
the  given  equipment  pair. 

Receiver  Response  Synthesis.  In  the  FAS  model,  the  receiver 
relative  response  characteristic  is  synthesized  by  four  line  seg- 
ments, which  are  linear  on  a logarithmic  scale,  and  appear  as 
shown  in  Figure  5.  The  variables  are  defined  below. 
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(18b) 


r(f)  = k , when  f <f  <f, 
s a 1 


(18c) 


1 or  f,  >f  >f_ 

D l 


r(f)  = k 


r(f)  = k 


where 


Then 


(v  - 

\W^l) 

(\(  - falY2 
\p^r/  • 


V when  f ^ < f 


(18d) 


when  f < f 


(18e) 


r(f)  = the  relative  response  of  the  receiver  at  frequency 

f. 


R(f)  = 


the  relative  response  in  dB  where  R(f)  * 10  log  r(f), 
and 

10  log  ks. 

lOn^,  10n2  respectively. 


Except  for  f ^ and  the  parameters  shown  are  required  inputs 

for  the  FAS  subroutine.  Methods  of  determining  these  parameters  are 
explained  below. 

Receiver  Tuned  Frequency,  fR.  The  tuned  frequency  of  the  receiver 

is  obtained  from  the  frequency  assignment  appropriate  to  the-  receiver 
being  considered,  and  is  a required  input'  to  the  FAS  program. 
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IF  Bandwidth,  B . The  intermediate  frequency  bandwidth  is 

obtained  from  the  nominal  characteristics  of  the  receiver,  usually 
set  forth  in  the  manufacturer's  data  in  technical  manuals  describing 
the  equipment,  or  it  is  obtained  from  measured  data.  It  is  a required 
input  and  must  have  a value  greater  than  zero. 

IF  Skirt  Slope,  N^.  This  parameter  is  extracted  from  the  given 

IF  selectivity  characteristics.  It  is  the  slope  of  the  skirt  in 
dB  per  decade.  For  example,  the  nominal  characteristics  of  a receiver 
may  specify  that  the  IF  selectivity  of  a receiver  has  a 20-dB  band- 
width of  BW  and  a 60-dB  bandwidth  of  BW2. 


Then  n^  is  obtained  from: 


(60-20)dB 

and 


40  dB  = lOrij 


log 


(19a) 


This  is  a required  input  and  must  be  greater  than  zero. 

RF  Skirt  Slope,  Characteristics  of  the  RF  circuitry,  which 

must  usually  be  theoretically  synthesized,  can  be  found  in  AVBASE, 
the  AVPAK  file  of  selected  avionics  equipment.  The  determination 
of  these  characteristics  is  the  result  of  special  analysis  intended 
to  study  such  characteristics.  There  are  two  methods  for  synthesizing 
RF  selectivity. 


The  first  method  requires  information  concerning  the  receiver 
design,  i.e.,  the  number  of  tuned  circuits  preceding  the  first  mixer 
in  a receiver.  This  information  can  be  obtained  by  examination  of 
the  circuit  diagrams  included  in  the  technical  manual  describing 
the  equipment.  When  the  number  of  tuned  circuits  is  known,  the 
relative  response  of  the  circuity  can  be  calculated  from  the  following 
equation: 
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= -lOn  log 


[■  -(^01 


(20b) 


where 


R p = the  RF  response  in  dB  relative  to  the  response  at 
the  tuned  frequency  of  the  receiver 

n = the  number  of  tuned  circuits  preceding  the  first 


f = the  frequency  at  which  the  relative  response  is 
' required 

Fr  = the  tuned  frequency  of  the  receiver 

) = the  selectivity  factor  of  each  tuned  stage,  defined 

S by 


where 


Rl  3 dB 


2 | f - f R | 3 b = t*ie  3 bandwidth  of  the  stage. 

The  parameter  Qs  is  not  generally  available  and  must  be  estimated. 
An  estimated  Qg  of  50  yields  sufficiently  conservative  results  for 
most  analyses. 

The  solution  of  Equation  20  as  a function  of  frequency  will 
yield  the  relative  selectivity  of  the  RF  circuitry  of  the  receiver. 

In  certain  cases,  the  number  ofj tuned  circuits  preceding  the 
mixer  is  unknown  but  the  RF  3-dB  bandwidth  is  specified  along  with 
the  image  rejection.  When  these  parameters  are  given,  an  alternate 
method  for  estimating  the  RF  characteristics  can  be  used. 

The  parameter  N2  is  equivalent  to  10n2,  which  equals  20n, 
where  n is  the  number  of  RF  tuned  circuits,  or  N2  = 10n2  = 20  n. 
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The  image  frequency  of  a receiver  is  separated  from  the  tuned 
frequency  by  twice  the  intermediate  frequency.  If  an  "image  band- 
width" is  defined  as  equal  to  four  times  the  intermediate  frequency, 
i.e.,  twice  the  image  frequency  separation  from  the  tuned  frequency, 
then  the  parameter  can  be  approximated  by: 


- 3 dB  = 10  n log 


fe)  or' 


(22a) 


- 3 dB 

TO 

V3  d J 


(22  b) 


where 

N2a  = the  ilPProximat:e  value  of  N2 
= the  image  rejection  in  dB 

BWj  = the  "image  bandwidth"  defined  above  (note  that  this 
bandwidth  is  a mathematical  device  and  not  a physical 
reality) 

BWj  = the  specified  3-dB  RF  bandwidth  of  the  receiver. 

The  model  does  not  accommodate  receivers  with  variable-tuned  inter- 
mediate frequencies. 

It  should  be  noted  that  the  universal  resonance  curve,  when 
plotted  on  a logarithmic  scale,  is  rounded  rather  than  linear  in 

the  vicinity  of  the  3-dB  bandwidth.  Therefore,  the  value  of  N_ 

# 2a 

obtained  from  Equation  22b  will  be  slightly  larger  than  is  appro- 
priate. Therefore,  this  value  should  be  rounded  off  to  the  nearest 

multiple  of  20  which  is  less  than  that  value.  Thus,  N_  = N„ 

*-  2a 

rounded  down  to  the  nearest  multiple  of  20.  When  this  value  of  N2 

is  determined,  then  n,  the  number  of  stages,  is  obtained  by  dividing 
by  20.  A corresponding  value  of  Q for  the  RF  circuitry  can  then 
be  calculated  by:  s 


(21'"  - .) 
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the  effective  selectivity  factor  for  each  tuned 
circuit  preceding  the  mixer 

the  overall  selectivity  factor  for  the  n tuned 
circuits 


n = the  number  of  tuned  circuits  as  determined  above. 

The  computation  of  the  two  parameters,  Q and  n,  results  in 

s 

values  which  can  be  used  in  Equation  20  to  estimate  the  RF  selec- 
tivity characteristics  of  the  receiver. 


If  neither  of  these  alternatives  is  feasible  due  to  a lack  of 
information,  the  preprocessor  will  assign  a worst-case  (i.e.  inter- 
ference susceptible)  recovery  value  of  20  dB/decade. 


Spurious  Response  Rejection,  The  minimum  spurious  response 

rejection  is  usually  specified  in  the  nominal  characteristics  of  the 
receiver.  If  not,  a worst-case  rejection  level  can  be  estimated  in 
the  following  manner.  Spurious  responses  arise  in  a superheterodyne 
receiver  when  a high  level  interfering  signal  combines  in  the  mixer 
circuitry  and  a product  at  or  near  the  intermediate  frequency  of 
the  receiver  is  generated.  In  general,  the  most  sensitive  of  these 
responses  arises  due  to  the  mixing  of  the  incoming  signal  and  the 
first  local  oscillator  frequency  in  the  first  mixer  stage.  The 
frequencies  at  which  the  interfering  signal  can  excite  these  responses 
is  given  by: 


Pflo  1 fIF 

q 


(24) 


where 

f = the  frequency  of  the  incoming  interfering  signal 

fj0  ■ the  injection  frequency  of  the  local  oscillator 

I 

f jp  ■ the  intermediate  frequency  of  the  receiver 
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The  fact  that  the  sign  preceding  fJF  in  Equation  24  can  take 

on  either  sense  (±)  indicates  that  for  a given  p,q  combination, 
a pair  of  responses  can  be  predicted,  one  for  the  negative  sense 
and  one  for  the  positive  sense.  In  a superheterodyne  receiver, 
the  local  oscillator  frequency,  flQ,  is  related  to  the  tuned 

frequency,  fR,  as  follows: 

flo  = fR  4 fIF  (25) 

where  the  positive  sense  is  appropriate  when  the  oscillator  fre- 
quency is  above  the  tuned  frequency  and  the  negative  sense  is 
appropriate  when  the  reverse  situation  occurs. 

The  combination  of  Equations  24  and  25  yields  two  sets  of 
relationships: 


when  flQ  > fR  (26a) 


when  fR  > flQ  (26b) 


These  expressions  enable  a determination  of  the  frequencies  at 
which  an  incoming  signal  can  result  in  the  most  sensitive  spurious 
responses.  Note  that  when  p»q=l,  the  two  responses  which  result 
are  the  receiver  response  to  its  tuned  frequency  and  the  response 
to  its  image  frequency.  The  relative  rejection  at  the  tuned  fre- 
quency is  zero  and,  from  the  previous  discussion,  it  is  known  that 
the  relative  rejection  at  the  image  frequency  is  merely  the  RF 
rejection  at  that  frequency. 


and 


sp 


sp 


(P+D  f 


p£r 

q 

and 


pfn  (P-1)  f 


IF 


IF 


pf^  _ (P+D  f1F 

q ’ q 

and 

Pf R (P-1)  fIF 
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The  relative  spurious  rejection  for  any  other  p,q  combination 

is  the  product  of  the  RF  rejection  at  the  incoming  frequency,  f , 

SP 

and  the  relative  mixer  conversion  loss  for  the  p,q  combination  being 
studied.  This  can  be  expressed  logarithmically  as: 

Ks  - Rr,,  ♦ uc  (27) 

where 

= the  relative  spurious  rejection  in  dB 

RRP  = the  relative  rejection  of  the  RF  circuits  at  f 
in  dB 

p = the  mixer  conversion  loss  of  the  actual  p,q  combina- 
tion relative  to  the  mixer  conversion  loss  of  the 
> p=q=l  combination. 

In  a previous  ECAC  measurement  effort  the  relative  mixer  conversion 
losses  for  transistor  mixers  were  established.  These  values,  which 
are  representative  of  most  mixers,  are  shown  in  TABLE  1. 


TABLE  1 

RELATIVE  TRANSISTOR  MIXER  CONVERSION  LOSS  IN  dB 


X 

1 

2 

3 

4 

5 

1 

0 

-65 

-76 

-84 

-83 

2 

-13 

-58 

-77 

-83 

-83 

3 

-13 

-65 

-77 

-81 

-82 

4 

-20 

-62 

-81 

-81 

-82 

5 

-22 

-62 

-78 

-84 

-82 

The  relative  spurious  rejection  level  obtained  with  Equation 
27  for  the  most  susceptible  response,  excluding  the  image  response, 
is  the  value  which  should  be  used  as  an  input  to  the  program.  If 
an  input  is  not  given,  a nominal  value  for  Kg  of  60  dB  is  assigned 

by  the  preprocessor.  The  image  rejection  will  also  be  essignad  a 
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nominal  value  of  60  dB  if  no  value  is  given,  but  the  image  response 
is  treated  as  a special  case.  When  this  situation  does  arise,  the 
receiver  is  synthesized  around  the  image  frequency  in  an  identical 
manner  to  the  synthesis  about  the  tuned  frequency,  except  that  the 
relative  threshold  is  reduced  by  the  input  image  rejection. 

Spurious  Response  Limit  Frequencies , f , f . As  stated  pre- 

A D 

viously,  these  frequencies  are  determined  from  the  point  where  the 
*CS  level  intersects  the  RF  selectivity  curve.  It  should  be  noted, 

however,  that  f^  and  f^  are  discrete  frequencies,  rather  than  fre- 
quency separations.  Since  the  RF  selectivity  is  usually  specified 
in  terms  of  response  versus  frequency  separation,  the  intersection 
point  is  more  readily  obtained  in  terms  of  a frequency  separation. 

It  is  necessary,  therefore,  to  add  this  separation  to  (or  subtract 
it  from)  the  receiver  tuned  frequency  to  obtain  the  appropriate 
values  of  ffl  and  f^.  If  discrete  frequencies  for  f#  and  f^  are 

not  entered  into  the  program,  recovery  values  will  be  assigned  by 
the  preprocessor.  For  the  fixed-tuned  case,  upper  and  lower  recovery 
values  are  determined  by  f ± 2IF.  For  the  range-tuned  case,  the 

lower  value  is  the  lower  limit  of  the  range  less  2IF  and  the  upper 
limit  is  the  upper  limit  plus  2IF. 

Transmitter  Spectral  Emission  Synthesis 

The  envelope  of  the  spectral  characteristics  of  a transmitter 
is  synthesized  by  three  line  segments  which  are  linear  on  a log- 
arithmic scale  as  shown  in  Figure  6. 


Figure  6.  Relative  spectral  emission  envelope  characteristic. 
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the  relative  level  of  the  spectral  emission  envelope 
at  frequency  f;  and, 

10  log  t(f),  the  relative  level  in  decibels. 


Each  of  the  transmitter  parameters  is  discussed  below. 

Conminications  Transmitters.  The  inputs  to  enable  synthesis 
of  the  spectral  envelope  of  communications  transmitters  are: 


B 


IT 


M, 


2T 


t, 


The  emission  envelope  bandwidth  of  the  first  break- 
points is  the  nominal  3-dB  bandwidth.  This  item 
has  a missing  input  recovery  mode  made  up  of  the 
FCC  emission  designator  and  pulse  characteristics 
derived  from  the  pulse  compression  indicator,  pulse- 
width,  and  pulse  rise  and  fall  times,  which  are 
all  required  inputs.  The  FCC  Emission  Designator 
is  made  up  of  the  bandwidth  containing  99%  of  the 
mean  radiated  power  and  the  modulation  type  used. 

The  first  slope  fall  off  of  the  emission  envelope; 
it  is  the  slope  adjacent  to  the  3-dB  bandwidth. 

The  value  given  must  be  >_  20  dB/decade.  This 
item  has  a missing  input  recovery  mode  similar  to 

B1T' 

This  is  the  emission  envelope  bandwidth  of  the 
second  breakpoints;  it  is  the  bandwidth  at  which 
the  envelope  shows  a second  fall-off  characteristic. 
This  item  has  a missing  input  recovery  mode  similar 

to  Bir 

M2  * This  is  the  second  slope  fall  off  of  the  emission 
envelope;  it  is  the  slope  at  frequencies  greatly 
separated  from  the  tuned  frequency.  The  value  used 
must  be  >_  40  dB/decade.  This  item  has  a missing 
input  recovery  mode  similar  to  B^. 

The  specification  of  the  minimum  20  dB/decade  for  M2  represents 

the  expected  minimum  fall-off  characteristic  of  the  noise  sidebands 
of  these  transmitters.  If  an  external  RF  filter  is  used  in  con- 
junction with  the  transmitter,  the  fall-off  characteristic  of  the 
filter  should  be  added  to  Mj. 
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Pulsed  Transmitters.  Inputs  for  pulsed  transmitters  can  be 
categorized  into  two  cases,  pulsed  transmitters  without  frequency 
modulation  during  the  pulse  interval  (PO  emission),  and  pulsed 
transmitters  with  frequency  modulation  during  the  pulse  interval, 
i.e.  "chirped"  pulsed  (P9  emission).  If  the  pulsed  transmitter 
concerned  has  a PO  emission,  the  pulse  width  (t)  entered  should 
be  the  width  at  the  half-amplitude  points.  If  the  emission  is 
P9,  the  v entered  should  be  the  "stretched"  pulse  width  or  total 
duration  of  the  pulse. 

For  PO  emissions,  the  spectral  envelope  can  be  found  using 
the  methods  described  by  Mason  and  Zimmerman. 14  The  results  are- 

t 

B1T  = 1 . 28/  (2t  ♦ tr  ♦ tf)  (29a) 


Mj  _>  20  dB/decade 


M2  40  dB/decade 


where 

T = fhe  pulse  width  between  one-half  amplitude  points 

ty  = the  time  required  for  the  pulse  to  rise  from  its 

10  percent  amplitude  point  to  its  90  percent  amplitude 
point,  i.e.,  the  pulse  rise  time 

J 

tf  * the  time  required  for  the  pulse  to  drop  from  its 

90  percent  level  to  its  10  percent  level,  i.e.,  the 
pulse  fall  time. 

4 

The  parameter  t is  a required  input  and  can  be  obtained  from 
the  nominal  characteristics  of  the  equipment.  The  parameters  B^, 

B2T*  “1.  tr*  tf  **®ve  *issing-input  recovery  modes. 

The  determination  of  the  envelope  parameters  of  a chirped 
pulse  can  be  quite  complex15  and  the  resulting  emission  spectral 


14Mason,  S.  and  Zimmerman,  H.,  Electronic  Circuits,  Signals,  and 
Systems,  John  Wiley  and  Sons,  New  York,  1960. 

lsKlauder,  J.  R.  et  al..  The  Theory  and  Design  of  Chirp  Radara, 
The  Bell  System  Technical  Journal,  July  1960. 
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envelope  for  such  a pulse  is  not  always  amenable  to  a three- line- 
segment  synthesis  technique.  However,  for  most  practical  eases 
the  inputs  can  be  determined  as: 


where 


(30a) 


tj  = the  total  frequency  deviation  during  the  pulse 

I)  = the  dispersion  ratio  = x f^  where  x is  the  total 
duration  of  the  pulse. 


(30b) 


Mj  = the  slope  of  the  line  in  dB/decade  which  joins  the 
®1T  ®2T 

points  ~y~  and  -j-  when  those  points  are  plotted 
on  a logarithmic  scale 

= second  slope  fall-off  of  the  emission  envelope. 

The  value  used  must  be  _>  20  dB/decade. 


The  inputs  Bjj»  ®2T*  ^1'  an<*  ^2  ^ave  roissing- input  recovery  modes. 


Harmonic  Emissions.  In  addition  to  synthesizing  the  envelope 
of  the  spectral  characteristicsof  the  fundamental  emission  of  the 
transmitter,  the  model  automatically  synthesizes  the  spectral 
envelopes  of  the  transmitted  harmonics  up  to  a maximum  of  the  ninth 
harmonic.  In  doing  so,  the  subroutine  assumes  that  the  spectral 
characteristics  of  each  harmonic  are  identical  to  that  of  the 
fundamental  except  that  the  reference  point  is  reduced  by  a level 
in  dB  equal  to  the  attenuation  specified  as  appropriate  to  the  har- 
monic. Attenuation  levels  for  all  harmonics  to  be  considered  (up 
to  the  ninth)  are  required  input  parameters. 

Rejection  Calculation 


After  the  receiver-response  and  spectral -emission  characteristics 
have  been  synthesized,  and  each  has  been  normalized  to  the  level 
appropriate  at  the  tuned  frequency  of  the  equipment,  the  expected 
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rejection  offered  by  the  receiver  to  the  undesired  emission  can 
be  calculated  using  only  frequency  relationships.  The  development 
of  this  concept  is  given  below. 

The  spectral  power  density  of  an  emission  can  be  approximated 


where 

PD  = the  spectral  density  in  watts/Hz 

p^  = the  average  transmitter  power  in  watts 

B = the  3-dB  emission  bandwidth  in  Hz. 

IT 

Thus  PD  represents  an  approximation  to  the  average  energy  con- 
tent in  the  emission.  However,  the  receiver  can  only  intercept  that 
energy  for  a period  of  time  equal  to  its  "response  time".  The 
"response  time",  tr,  of  a receiver  can  be  considered  to  be  the  inverse 

of  its  3-dB  bandwidth.  Accordingly,  the  maximum  power  transfer  between 
a transmitter  and  receiver  can  be  calculated  by: 


r tR  bit  (tr} 

where 


P = the  average  received  power  in  watts 

B = the  3-dB  bandwidth  of  the  receiver, 

r 

The  receiver  rejection,  lf,  is  the  ratio  of  the  received  power  to 
the  transmitted  power: 

P Bd 

1 - _E.  s —5—  or  in  logarithmic  terms,  (33) 

f pt  b1t 


L-  ■ 10  log  1-  = 10  log  „ — 

£ f B1T 
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The  maximum  power  transfer  occurs  when  the  transmitter  and  receiver 
are  tuned  to  the  same  frequency.  Therefore,  a cochannel  rejection 
factor  (lf)CQ  can  be  defined  as: 

<Vco  * 17;  "hen  “it  - br  <54> 

and 

1,  when  B1t  <_  BR 


The  reason  for  the  two  cases  is  evident  when  it  is  remembered 
that  a receiver  cannot  receive  more  power  than  is  transmitted. 

It  should  be  noted,  however,  that  the  cochannel  rejection 
factor  in  Equation  34  has  been  defined  in  terms  of  average  power. 
Wien  the  undesired  emission  is  from  a pulsed  transmitter,  the  peak 
received  power  is  usually  of  more  interest  than  the  average  received 
power.  Thus  for  a pulsed  transmitter: 


‘V.vg 

a 

(PT)pk  t (PRF) 

where 

‘Vavg 

s 

the  average  power  in  watts 

(Vpk 

s 

the  peak  power  in  watts 

T 

s 

the  pulse  duration  in  seconds 

PRF 

s 

the  pulse  repetition  frequency  in  hertz. 

If  Equation 

35 

is  combined  with  Equation  32  then: 

<Vavg 

s 

(ptV  ’ <™F>  577 

(35) 


(36) 


However,  in  a pulsed  transmitter,  t * , and  the  average  received 

B1T 

power  (pR)avg  is  related  to  the  peak  received  power,  (Pr^*  by: 
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where 

tr  = the  response  time  of  the  receiver,  which  is  approx- 
imately pi-  . 


Thus,  combining  Equation  37  with  Equation  36  and  defining  a cochannel 
rejection  factor  for  pulsed  emissions,  it  is  seen  that: 


I 1 , when  B1T  <_  By 

Equation  38  applies  for  considerations  involving  the  peak  power 
transfer  due  to  a pulsed  transmitter. 


Since  the  cochannel  rejection  has  been  determined,  the  total 
rejection,  Lf,  at  any  frequency  can  be  calculated  by: 


Lf  " (Lf)co  + R(fT}*  W 

Equation  39  applies  for  considerations  involving  the  power  transfer 
resulting  from  the  fundamental  emission  of  the  transmitter  as 
affected  by  receiver  selectivity. 

Also, 

Lf  " (Lf>co  * T(fR)  («0) 

applies  for  considerations  involving  the  power  transfer  resulting 
from  emission  sidebands  which  occur  within  the  receiver  passband. 

In  Equations  39  and  40, 

R(fT)  ■ the  relative  response  of  the  receiver  at  frequency 

fT 

T(fR)  ■ the  relative  emission  level  of  the  transmitter  at 
frequency  fR. 
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The  FAS  subroutine  essentially  calculates  the  value  of  Lf  for 

each  of  the  situations  given  above,  compares  the  two  values,  and 
selects  the  lowest  of  the  values  obtained  for  use  in  Equation  7. 
However,  since  both  the  relative  response  characteristic  of  the 
receiver  and  the  relative  emission  level  of  the  transmitter  can 
vary  over  a range  of  frequencies,  the  computation  is  made  using 
integration  techniques.  ■, 

Further  details  on  the  subroutine,  including  a rigorous  math- 
ematical description  of  the  calculation  techniques,  may  be  found 
in  Reference  13. 

Degradation  Thresholds,  (S/I)^. 


The  input  signal-to-interference  ratio  at  which  operational 
degradation  begins  to  occur  in  a receiver  is  defined  as  the  degra- 
dation threshold,  and  is  identified  by  the  symbol  (S/I)T.  It  is 

the  minimum  signal-to-interference  ratio  corresponding  to  non- 
interference and  may  be  obtained  from  equipment  characteristics. 

If  no  ratio  is  given,  a value  of  20  dB  is  assigned  by  the  pre- 
processor. 

POWER  DENSITY  CALCULATION 


The  capability  to  predict  power  densities  resulting  from 
transmitters  located  on  an  airframe  as  well  as  from  transmitters 
located  on  neighboring  aircraft  is  one  of  the  features  AVPAK  5 
provides. 

The  power  density  at  a given  point  resulting  from  the  operation 
of  a single  transmitter  is  calculated  by: 

PD  = PT  + Gt  - Lp  ♦ 20  log  f - 38.5  (41) 

where 

PD  = power  density  level  in  dBm  per  square  meter 
PT  * transmitter  output  power  in  dBm 

Gj  * effective  transmitter  antenna  gain  in  the  direction 
of  the  point  in  dBm 

L ■ coupling  path  loss  between  transmitting  antenna  and 
p desired  point  in  dB 

f * the  transmitting  frequency  in  MHz. 
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This  equation  may  vary  for  different  transmitter  modulation 
types.  AVPAK  3 considers  two  different  modulation  types,  pulsed 
and  non-pulsed  (continuous  wave).  For  pulsed  transmitters,  (mod- 
ulation type  P9  or  PO) , both  peak  and  average  power  densities 
are  calculated.  For  non-pulsed  transmitters  an  average  power 
density  is  used.  The  revised  power  density  expressions  are: 


p 

rDPK 

= 

PTPK  + GT  - Lp  + 20  lo8  f - 38.5 

(42) 

PDAVG 

= 

PTAVG  * ^T  ' kp  + 2G  1°8  ^ ” 38.5 

(43) 

where 

p 

DPK 

= 

the  peak  power  density  level  in  dBm  per  square 
meter 

PDAVG 

= 

the  average  power  density  level  in  dBm  per 
meter 

square 

p 

TPK 

s 

peak  transmitter  output  power  in  dBm 

PTAVG 

= 

average  transmitter  output  power  in  dBm. 

The  transmitter  output  power  for  non-pulsed  transmitters  may 
be  obtained  from  nominal  characteristics  or  measured  data.  This 
value  is  assumed  to  be  d.n  average  power  level.  For  pulsed  trans- 
mitters the  output  power,  obtained  in  the  same  way,  is  assumed  to 
be  a peak  power  level.  It  is  necessary  to  calculate  the  average 
power  output  level  for  pulsed  transmitters.  This  is  done  using 
the  following  formula: 


Ptavg  " (PtPk>  (T)  (PRF) 

where 


(44) 


tavg 

Ptpk 

x 


PRF 


average  transmitter  output  power  in  kilowatts 
peak  transmitter  output  power  in  kilowatts 
transmitter  pulse  width  in  seconds 
transmitter  pulse  repetition  frequency  in  pps. 


AVPAK  3 also  calculates  the  cumulative  power  density  due  to 
the  effects  of  all  transmitters  at  each  point  of  interest,  following 
the  individual  transmitter  calculations.  The  cumulative  power 


FAA-RD-76- 50 


Section  3 


density  uses  the  average  power  densities  so  that  both  pulsed  and 
non-pulsed  transmitters  are  represented.  It  is  calculated  as: 


DSC 


i=l 


(45) 


where 


DSC 


DAVG  i 


cumulative  average  power  density  due  to  all  trans- 
mitters, at  point  of  interest,  dBm  per  square  meter 

number  of  individual  transmitters  in  the  environment 

average  power  density  at  point  of  interest  due  only 
to  transmitter  i in  dBm/ square  meter. 


The  field  may  also  be  expressed  in  terms  of  the  magnitude  of 
the  electric  field  vector  in  volts  per  meter.  The  conversion  in 
the  program  is  done  thusly: 


= 10 


(46) 


where 

Fg  = field  strength  in  volts/meter 
PD  = power  density  in  dBm/ square  meter. 

Both  the  power  density  and  electric  field  strength  are  computed 
and  printed.  Note  that  field  strength  is  computed  for  the  case 
where  Pp  is  equal  to  Ppg^.  (cumulative  power) . The  value  so  computed 

is  fictitious  and  could  never  be  verified  by  direct  measurement. 

It  represents  an  equivalent  field  strength  which,  if  produced  by 
a single  transmitter,  would  have  the  same  effect  as  the  sum  of 
several  discrete  sources. 

A conversion  chart  from  dBm/m2  to  volts/meter  is  shown  in 
Figure  7. 
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strength  in  volts/meter. 
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Coordinate  Systems 

Coordinates  for  transmitters,  receivers  and  points  may  be 
entered  in  the  model  by  the  user  in  rectangular  (x,y,z)  coordinates, 
cylindrical  (p,  0,  z)  coordinates,  or  the  buttl ine-waterl ine- 
fuselage  station  (B,  W,  F)  coordinate  system. 

Regardless  of  what  coordinate  system  is  used,  the  location 
coordinates  are  converted  in  the  model  to  the  rectangular  system. 

Referring  to  Figure  S,  in  the  rectangular  coordinate  system 
the  origin  is  forward  at  the  center  of  the  fuselage  with  the  pos- 
itive X direction  being  to  the  right  (starboard)  side  of  the  air- 
craft. The  origin  is  in  a similar  location  for  the  cylindrical 
coordinate  system.  In  buttline-watcrline-fuselage  station  coordi- 
nates the  origin  is  forward,  not  at  the  center  but  at  the  bottom 
of  the  fuselage.  The  positive  buttline  is  to  the  left  (port) 
side  of  the  aircraft,  the  waterline  direction  toward  the  top  of 
the  aircraft,  and  the  fuselage  station  direction  toward  the  tail 
of  the  aircraft. 

OUTPUT  TYPES 


AVPAK  3 provides  three  ways  of  solving  the  basic  interference 
expression  to  determine  if  degradation  is  to  occur  between  a 
transmitter  and  a receiver.  The  three  analysis  types  are  deter- 
ministic, functional  and  probabilistic. 

Deterministic  Analysis 

A deterministic  analysis  is  performed  in  AVPAK  3 by  considering 
specific  transmitter/receiver  pairs  and  evaluating  the  following 
equation: 

Pld  = PT  ♦ GR  ♦ gt  - LP  ♦ (S/I)T  ‘ RS  (47) 

where  all  terms  have  been  previously  defined. 

If  Plcl  is  greater  than  zero,  then  AVPAK  3 makes  the  prediction 

that  interference  to  the  receiver  due  to  the  transmitter  will  occur, 
and  the  particular  equipment  pair  is  entered  in  the  list  of  potential 
interference  cases.  If  P^  is  less  than  zero,  the  equipment  pair 

is  entered  in  the  non-interference  list. 

The  values  of  the  terms  in  Equation  47  are  based  on  nominal 
(generally  median)  equipment  parameters.  The  user  may  either  enter 
these  equipment  parameters  or,  by  using  a unique  identification  code. 
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automatically  retrieve  the  desired  parameters  from  the  AVPAK 
file  AVBASE. 

Functional  Analysis 

A functional  analysis  can  be  performed  instead  of  a purely 
deterministic  one.  This  type  of  analysis  is  identical  to  the 
deterministic  type  except  a different  data  base,  AVFILE  is  called 
upon  for  equipment  parameter  values.  In  AVFILE,  equipments  are 
grouped  according  to  functional  classes.  A functional  class  is 
a group  of  equipments  whose  characteristics  are  similar  enough 
so  that  all  individual  equipments  in  the  group  can  be  represented 
bv  a single  representative  equipment  type.  The  parameter  values 
which  represent  a single  representative  equipment  type  are  the 
median  values  of  the  equipments  making  up  the  functional  class. 

Probabilistic  Analysis 

The  third  type  of  analysis,  probabilistic,  allows  the  trans- 
mitter/receiver pairs  to  be  analyzed  using  statistical  descriptions 
or  cumulative  distributions  for  the  terms  in  Equation  47  whose 
expected  errors  from  median  values  can  be  confidently  described. 

In  other  words,  median  values  for  P^.  (transmitter  power)  Lp  (off- 

frequency  rejection)  Lp  (path  loss)  and  Rg  (receiver  sensitivity) 

may  not  be  realistic  enough  to  yield  an  accurate  prediction  of 
interference  using  a deterministic  approach.  The  distributions 
for  these  parameters  are  stored  along  with  the  functional  class 
data  in  AVFILE. 


The  median  values  for  the  terms  in  Equation  47  are  denoted 
as  PT,  Lp,  Lp,  and  Rg.  The  expected  variation  of  these  median 
values  are  epT,  eT0FR  and  eR0FR,  eLp,  and  eRS>  eTOFR  and  eRQFR 

represent  the  variation  in  transmitter  and  receiver  off-frequency 
rejection,  respectively. 

The  probabilistic  analysis  begins  with  a procedure  similar 
to  the  deterministic  analysis.  The  following  inequality  is  checked: 


P,N  * DU,  > Rs 


where 


(48a) 


PIN  “ PT  + GT  * GR  * LP  - 4 * 


(48b) 
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and 


f 


IN 


the  effective  median  interference  power  for  the 
functional  classes  being  analyzed  (dBm). 


R_  = median  receiver  sensitivity 


DU, 


■ [<“W2  * ‘“"LP)2  * (0UT0FR)2  * '““rOPr)2 


where 


<DlW2 


T 


(48c) 


DU, 


upper  decile  (90%  probability)  correction  factor 
in  dB  made  up  of  upper  decile  values  determined 


from  cumulative  error  distributions  of  P, 


up. 


Lp  and  Rg.  The  Lp  factor  is  comprised  of  the 

contribution  of  the  transmitter  to  the  receiver 
off- frequency  rejection  (TOFR)  cumulative  error 
and  the  contribution  of  the  receiver  to  the 
receiver  rejection  cumulative  error  (ROFR). 

These  upper  decile  values  are  determined  from  the  cumulative  error 
distribution  plots  of  the  respective  equipment  parameters  at  the 
90%  probability  level. 

If  Expression  48a  is  not  satisfied  by  an  equipment  pair,  the 
pair  is  simply  entered  in  the  list  of  predicted  non-interference 
cases. 


If  Expression  48a  is  found  to  be  true  for  a given  pair,  the 
probabilistic  analysis  continues  by  statistically  convolving  the 
cumulative  distribution  of  expected  error,  for  each  of  the  five 
factors.  The  convolved  distribution  for  any  particular  equipment 
pair  is  determined  as: 


:PT  * eT0FR 


'ROFR 


'LP 


'RS 


(49) 


where 


i 

eT 


convolved  cumulative  error  distribution  function 
for  the  transmitter/ receiver  pair,  expressed  as 
dB  error 
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e 


PT 


CTOFR 


eROFR 


c 


LP 


convolved  cumulative  distribution  of  expected 
errors  from  the  median  value  of  transmitter  power 
in  dB 

convolved  cumulative  distribution  of  expected 
errors  from  the  median  value  of  the  transmitter 
contribution  to  the  receiver  off-frequency  rejec- 
tion, in  dB 

convolved  cumulative  distribution  of  expected 
errors  from  the  median  value  of  the  receiver 
contribution  to  the  receiver  off-frequency  rejec- 
tion, in  dB 

convolved  cumulative  distribution  of  expected 
errors  from  the  median  value  of  path  loss,  in  dB 


* convolved  cumulative  distribution  of  expected 

errors  from  the  median  value  of  receiver  selectivity 
in  dB. 


The  terms  TOFR  (transmitter  contribution  to  receiver  off-fre- 
quency rejection)  and  ROFR  (receiver  contribution  to  receiver  off- 
frequency  rejection)  require  additional  development.  Their  values 
and  associated  deviations  are  obtained  in  the  following  manner. 

The  distributions  of  TOFR  and  ROFR  are  assumed  to  be  normal,  and 
the  median  values  for  both  terms  are  set  equal  to  zero.  The  stand- 


ard deviation  of  ROFR,  a 


ROFR’ 


is  associated  with  the  discrete 


spurious-response  level,  Kg,  discussed  previously  under  Receiver 

Rejection,  page  31.  When  groupings  of  individual  equipments  were 
taken  to  form  the  functional  classes,  the  median  value  and  standard 

The 


deviation  of  Kg  was  determined  for  each  functional  class. 


value  of  a 


ROFR 


was  set  equal  to  the  associated  Kg  value. 


A cumula- 


tive distribution  of  expected  errors  from  the  median  value  was  then 


formed  and  convolved,  yielding  t 


ROFR’ 


The  standard  deviation  for 


TOFR,  °TQpR»  was  determined  through  transmitter  emission-envelope 

synthesis.  For  each  transmitter  functional  class,  the  spectral 
emission  envelope  was  modeled  using  median  values.  The  standard 
deviations  of  the  parameters  used  were  known.  The  value  of  the 
appropriate  o was  then  added  to  the  median  value  of  the  emission 
bandwidth  defined  on  page  39)  at  which  the  envelope  fall- 


off  slope  changes  (i.e.,  B-_  ♦ oD  ).  The  vertical  difference 

2T  B2t 
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between  the  envelope  modeled  using  the  median  value  of  B2T  and  the 
envelope  modeled  using  (B^  ♦ 0g  ) , at  the  beginning  of  the  second 


slope,  is  one-half  the  value  of  a, 


TOFR" 


It  is  one-half  the  value 


because  the  total  variation  about  B2T  would  also  contain  the  area 

one  sigma  below  B^  (i.e.,  B2T  - o_  );  because  of  symmetry  it  is 

2T 

only  necessary  to  calculate  one  case.  After  values  of  the  TOFR 
parameter  have  .been  obtained,  a cumulative  distribution  is  formed 
and  convolved  to  arrive  at  eT0FR' 

The  convolved  cumulative  error  distribution,  eT>  is  then  plotted 

in  the  computer  output  as  the  predicted  probability  of  interference 
P(I)  versus  dB  error  (e).  Since  the  input  distributions  are  repre- 
sented by  worst-case  interference  conditions  at  the  high  end  of 
the  probability  scale  for  the  most  positive  cumulative  error  values, 
the  convolved  distribution  then  represents  the  worst-case  value 
in  a similar  manner.  The  distribution  type  may  be  normal,  or  any 
other  type,  which  may  be  determined  by  inspecting  the  plot.  The 
probability  of  interference  for  a given  pair  may  then  be  determined 
by: 


e ■ PIN  ' RS 


where 


* abscissa  coordinate  value  of  the  plot  eT,  from 

which  the  desired  ordinate  P(I),  (predicted  prob- 
ability of  interference)  can  be  determined. 


As  an  example,  see  Figure  9.  Using  Equation  50,  if  is 

IN 

determined  to  be  -95  dBm  and  R has  a value  of  -100  dB,  c is  equal 

O 

to  5 dB.  From  inspection  of  the  plot  in  Figure  9 one  can  deter- 
mine the  probability  of  interference  for  this  transmitter/receiver 
pair;  there  is  a 70%  probability  of  interference.  It  can  als»  be 
seen  from  the  graph  that,  to  decrease  the  probability  of  inter- 
ference to,  for  example,  50%,  the  value  of  PjN  must  be  decreased 
5 dB  to  -100  dB. 


52 


FAA-RD-76-S0 


Section  3 


DATA  BASES 

There  are  two  data  base  files  associated  with  AVPAK  3: 

1.  AVRASE,  which  carries  equipment  data  in  a format 
necessary  to  execute  AVPAK  in  the  deterministic  mode,  and 

2.  AVFILE,  which  carries  functional  characteristics 
necessary  to  execute  AVPAK  in  the  probabilistic  mode. 

In  addition,  all  equipment  data  collected  in  the  course  of  these 
projects  has  been  placed  in  the  ECAC  Nominal  Characteristics  File 
(NCF) . 

A description  of  AVBASE  and  AVFILE  is  presented  below.  Both 
files  have  been  augmented  with  new  avionics  equipment  records 
which  have  been  obtained  since  AVPAK  2 was  completed.  The  contents 
of  AVBASE  and  AVFILE  appear  in  APPENDIX  D. 

AVFILE 

The  functional  class  probabilistic  file,  AVFILE,  stores  equip- 
ment characteristics  for  distinct  functional  classes  of  equipment, 
as  defined  earlier  in  this  section.  There  are  four  functional 
transmitter  classes;  VHF  communications,  ATCRBS  transponder.  Weather 
Radar,  and  DME/TACAN.  Altimeters  have  been  removed  as  a functional 
class  because  of  limited  data  available.  They  are  retained  in 
AVBASE.  There  are  six  functional  receiver  classes:  VHF  communi- 

cations (of  which  there  are  three  subclasses),  VOR  (VHF  omni-range)/ 
Localizer  (of  which  there  are  also  three  subclasses) , Glideslope, 
ATCRBS  transponder,  Weather  Radar,  and  DME/TACAN.  For  both  the 
VHF  communications  and  the  VOR/Localizer  classes,  subclass  1 con- 
tains those  equipments  with  maximum  spurious  responses  that  are 
as  much  as  89  dB  down  from  the  on- tune  response,  subclass  2 from 
90-109  dB  down,  and  subclass  3 from  110-140  dB  down.  These  two 
functional  classes  have  subclasses  because  of  their  large  standard 
deviations;  by  subclassifying  them  the  standard  deviations  have 
been  significantly  reduced. 

Each  class,  represented  by  a unique  one-digit  or  two-digit 
functional  code,  contains  basic  parameter  data  and  a description 
of  the  sampled  data  used  to  derive  the  parameters.  Numeric  codes 
1 through  49  inclusive  are  reserved  for  transmitters  and  codes 
50  through  99  inclusive  are  reserved  for  receivers.  There  are 
no  functional  antenna  classes. 

AVBASE 

Execution  of  AVPAK  in  the  deterministic  mode  requires  para- 
meters not  readily  available  from  technical  manuals.  Specific 
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analyses  of  technical  manuals  was  accomplished  at  ECAC  in  order 
to  develop  such  characteristics  as  IF-  and  RF-slope  fall  off s,  etc. 
AVBASE  may  be  considered  as  a storage  vehicle  for  these  character- 
istics. It  should  be  noted  that  these  values  have  been  specifically 
developed  to  drive  the  AVPAK  model  and  are  not  for  general  use. 

Those  values  found  in  the  ECAC  Nominal  Characteristics  File  are 
for  general  usage. 

Each  AVBASE  record  is  represented  by  a unique  identifier 
code  (see  TABLE  2) . Numeric  codes  5100  through  8999  inclusive 
are  reserved  for  transmitters,  codes  100  through  3999  inclusive 
are  reserved  for  receivers. 


TABLE  2 


AVBASE  EQUIPMENT  IDENTIFICATION  NUMBERS 


Manufacturer 

Receivers 

Transmitters 

ARC-CESSNA 

100 

5100 

BENDIX 

500 

5500 

COLLINS 

900 

5900 

DYNAIR  - RADAIR 

1300 

6300 

EDO  - AIRE 

1700 

6700 

GENAVE 

2100 

7100 

KING 

2500 

7500 

NARCO 

2900 

7900 

RCA 

3300 

8300 

WILCOX 

3700 

8700 

Function 

Receiver  Sub-I.D. 

Transmitter  Sub-I.D. 

VHF  Omni-range/ localizer 

00  - 49 

VHF  Communications 

50  - 99 

50  - 99 

Glide  Slope 

100  - 139 

ATCRBS  Transponder 

140  - 179 

Weather  Radar 

180  - 219 

IfMlfi 

Distance  Measuring  Equipment 

220  - 239 

Altimeter 

260  - 299 

Note:  An  equipment  I.D.  is  created  by  adding  the  manufacturer’s  I.D. 

to  the  functional  Sub-I.D.  For  example,  a Bendix  communications 
receiver  number  would  fall  between  550  and  599  inclusive. 
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PROGRAM  UTILIZATION 


GENERAL 


This  section  discusses  the  control  and  data  cards  that  are 
necessary  to  run  AVPAK  3.  APPENDIX  E contains  several  sample 
runs  for  the  various  calculations  and  configurations  that  are 
handled  by  the  program. 

BASIC  AVPAK  3 DECK 


The  standard  card  deck  used  in  executing  the  AVPAK  3 model 
on  ECAC's  UNIVAC  1110  Computer  is  as  follows: 

@RUN 

0QUAL  HASELTINE 
@MAP,  IS,  AVPAK3 
IN  ECAC*  L IF I L/U . AVPAK3 
LIB  ECAC*LIFIL/U. 

NOT  ECAC  * L I F I L/U . RANTEN 
LIB  ECAC*MODLIB/U. 

LIB  ECAC*LIB/U. 

0SEC,  U 
@XQT  AVPAK3 

(Data  Sets) 

0EOF 

0PMD,E 

0FIN 

DATA  CARDS 


The  card  formats  for  the  data  cards  are  listed  in  TABLES  3 
through  12  and  the  coordinate  system  for  determining  main  beam 
pointing  angles  in  intersite  analyses  is  illustrated  in  Figure  10. 
Details  on  these  cards  and  also  on  program  utilization  can  be 
found  in  Reference  12. 
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TABLE  3 

GENERAL  PARAMETER  DATA  CARD 


Required 

Coluan(s) 

Format 

Description 

Units 

Yes 

1 - 3 

INT* 

Total  niaber  of  transmitters  (50  Maximum) 

Yes 

4-6 

I NT® 

Total  number  of  receivers  or  power  density  points  (50  maximum) 

Yes 

7 - 10 

INT* 

Total  number  of  antennas  (100  maximum) 

Yes 

12 

INT* 

Calculation  desired:  (0  or  8 * interference  to  noise  ratio, 

1 » power  density) 

Yes 

14 

INT* 

Analysis  desired:  (0  or  8 * cosite,  1 • intersite) 

Yes 

16 

INT* 

INR  answer  type:  (0  or  8 • not  applicable,  1 • deterministic, 

2 • functional,  3 ■ probabilistic) 

Yes 

18 

INT* 

Coordinate  system  to  be  used  for  inputs:  [1  ■ rectangular 

(X,  Y,  Z),  2 » buttline,  waterline  fuselage  station  (B,  N,  F), 
3 • cylindrical  (p,  8,  Z)] 

Yes 

19  - 24 

FPb 

Maximum  fuselage  radius 

in. 

No 

25  - 50 

FPb 

Bulkhead  distance  aft  from  nose 

in. 

No 

31  - 36 

FPb 

Bulkhead  height  above  fuselage  centroid 

No 

40  - 78 

FLDC 

Title  of  run 

aINT:  integer  value,  must  be  right  justified. 

^FP:  floating  point  value,  must  contain  a decimal  point. 

CFLD:  field  data  (character  data). 
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TABLE  9 

TRANSMITTER  DATA  CARDS 


Required 

Dl* 

D2b 

FC 

p" 

Column(a) 

Format 

Description 

Units 

Yes 

Yet 

Yet 

Yet 

1 - 3 

1NT 

Number  of  antennas  used  by  transmitter 

No 

No 

No 

No 

4 - IS 

FID 

Nomenclature 

Yet 

No 

No 

NO 

17  - 22 

FP 

lover  operating  frequency 

MHz 

Yet 

No 

No 

No 

24  - 29 

FI* 

Upper  operating  frequency 

MHz 

No 

Yet 

Yet 

Yet 

JO  - JJ 

(NT 

Functional  file  (AVFILE)  or  Oata  8a se 
file  (AVBASE)  code  number 

NO 

No 

No 

No 

JS  - 58 

FP 

First  cmissien  bandwidth 

Ml 

No 

No 

No 

No 

40  - 43 

FP 

Second  emission  bendvidth 

MHz 

No 

No 

No 

NO 

45  - 48 

PP 

First  emission  slope  falloff 

dB/decade 

No 

No 

No 

No 

SO  - S3 

FP 

Second  emission  slope  falloff 

di/decade 

Yet 

No 

NO 

No 

55  - 57 

FP 

Power 

dhn 

r« 

No 

NO 

No 

58  - 6 2 

FP 

Emission  designator  > bandwidth 

MHz 

Yet 

No 

NO 

No 

64  - 65 

FLO 

Emission  designator  - modulation  type 

No 

No 

No 

NO 

67 

FLU 

Pulse  compression  indicator  "BManot 
pulsed,  "c"-  pulsed 

Yet 

No 

No 

No 

68  - 71 

FP 

Pulse  width 

I4S 

Yet 

No 

No 

No 

72  - 75 

FP 

Average  pulse  rise  and  fall  time 

MS 

No* 

NO 

No 

No 

76  - 79 

FP 

Pulse  repetition  frequency 

kHz 

NO 

Yet 

To. 

Yet 

80 

INT 

Presence  of  second  transmitter  data 
card  indicator,  B-no  second  card 
present,  all  data  on  first  card,  1-read 
second  card  for  additional  transmitter 
data 

second 

card 

NO 

No 

No 

No 

1 - 6 

FP 

lower  transmitter  filter  limit 

Wz 

No 

No 

No 

NO 

8 - 13 

FP 

Upper  transmitter  filter  limit 

Yet 

NO 

No 

NO 

15 

INT 

Number  of  transmitter  harmonics  to 
consider  in  EMC  analysis  (1,  2,  3 9) 

NO 

No 

No 

NO 

17  - 20 

FP 

Suppression  level  of  second  harmonic 

dl 

No 

No 

No 

No 

22  - 25 

FP 

Suppression  level  of  third  harmonic 

dB 

No 

No 

No 

NO 

27  - JO 

FP 

Suppression  level  of  fourth  harmonic 

dB 

No 

NO 

No 

No 

32  - 55 

FP 

Suppression  level  of  fifth  harmonic 

dB 

NO 

No 

No 

No 

37  - 40 

PP 

Suppression  level  of  sixth  harmonic 

dB 

No 

No 

No 

No 

42  - 45 

FP 

Suppression  level  of  seventh  harmonic 

dB 

No 

NO 

No 

No 

47  - SO 

FP 

Suppression  level  of  eighth  harmonic 

dB 

No 

No 

No 

NO 

$2  - 55 

FP 

Suppression  level  of  ninth  harmonic 

dl 

aDl:  deterministic  type  run,  not  using  AVBASE. 

bD2:  deterministic  type  run,  using  AVBASE. 


F:  functional  type  run. 

^P:  probabilistic  type  run. 
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TABLE  10 


COMENTS  ON  TRANSMITTER  DATA  CARDS 


Field  or  P«r«Mt«r 

Acceptable  Range  or  Value 

Recovery  Value 

NMonclature 

N/A 

A 12  character  randou  dummy  tmmm 

PuIm  coapmtion  indicator 

Must  be  blank  or  "C" 

Set  to  "CM  if  Modulation  typo  it 
P9 

Pulse  rise/ fall  tiMO 

Must  be  less  than  pulsewidth 

1/10  of  pulsewidth,  PM 

First  Mission  bandwidth 

Must  be  less  than  second  enission 
bandwidth;  product  of  first  eeiission 
bandwidth  and  puliewidth  aust  be  at 
1.0  for  a chirped  equipment 

For  Modulation  typo  A-  or  F-  sot 
to  Mission  designator  bandwidth 
EDM.  For  Modulation  typo  F9,  sot 
to: 

(EDM)  ^.^T0aB|_rwy  ^ 

For  Modulation  typo  PO,  sat  to: 

1.28 
(2)”  ■(*) 

Second  eaission  bandwidth 

Must  be  greater  than  first  eaission 
bandwidth 

For  aodul.tioo  type  A-  or  F-  sat 
to  EDBN  tiaos  10.  For  oodoUtioo 
type  P9,  set  to: 

(EDM)  (1*Vi»w5  m) 

For  Modulation  typo  PO,  sot  to: 
(pulse  rise/fall  tiuo) 

First  eaission  slope  fallofl 

>20  dB/decade 

For  modulation  type  A*  or  F-  sot 
to  SO  di/docado.  For  aodolatioo 
typo  PO  sot  to  20  db/docado.  For 
oodulation  type  PD  aa  approslmottoo 
is  aad.  using  PN,  puis,  riso/foll 
tine,  and  firat  aad  eecoad  amiaaioa 
bandwidth. 

Second  eaission  slope  fall- 
off 

> 20  dt/de cade 

For  Modulation  types  A*  or  P- 
sot  to  20  dl/doondo.  For  Mod- 
ulation typos  PO  or  P9  sot  to 
40  dh/docodo. 

Makar  of  hasmooics  to  bo 

M and  <# 

Sat  to  1 if  blank 

Harmonic  suppression  levels 

N/A 

Sot  to  60  dl  if  blank 

(Relating  fro+iaary 

Upper  frequency  aust  be  greater  than 
or  equal  to  lower  frequency,  average 
operating  frequency  aust  be  at  least 
SO  Wit 

N/A 

Filter  fn^ii airy  limit 

Upper  frequency  oust  be  greater  than 
or  equal  to  lower  frequency 

N/A 

Pulsewtdth 

Miat  ba  g re. tar  than  pulaa  rita/fall 
time;  product  of  firat  eaission  bond’ 
width  and  puli  .width  aunt  bo  at  laaot 
1 for  a chlrpod  onuipaoat. 

N/A 

i i 
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SECTION  5 
SUMMARY 


GENERAL 


AVPAK  is  intended  as  a tool  for  identifying  antenna-coupled 
interference  problems  on  board  aircraft.  It  will  compute  inter- 
ference levels  in  avionics  receivers  from  transmitters  located 
on  the  same  aircraft  or  at  points  remote  from  the  aircraft. 
Further,  it  will  compute  power  densities  at  any  designated  point 
from  a transmitter  or  a number  of  transmitters. 

AVPAK  does  not  consider  near- field  antenna  anomalies,  nor 
does  the  model  consider  high-power  nonlinear  effects. 

MODEL  CAPABILITIES 

The  composite  capabilities  of  the  AVPAK  3 model  are  listed 
and  discussed  below,  including  responses  to  the  most  recent 
objectives  set  forth  in  Inter-Agency  Agreement  FA  70  WAI-175, 

Task  Assignment  No.  26. 

1.  Calculation  of  the  interference  effects  between 

both  cosite  and  intersite  transmitters  and  receivers  operating 
on  an  aircraft.  This  is  done  by  determining  the  expected  level 
of  interference  relative  to  the  degradation  threshold  of  each 
receiver.  If  this  expected  level  equals  or  exceeds  the  thresh- 
old value,  degradation  is  expected  to  occur.  This  calculation 
will  take  into  account  the  following  factors:  antenna  paths 

including  those  that  encounter  bulkheads,  airfoils,  and  wing 
weapon  pods,  and  also  between  antennas  which  are  raised  from  the 
fuselage  skin.  As  an  additional  user  convenience,  a preprocessor 
will  insert  certain  nominal  predetermined  values  if  certain  trans- 
mitter or  receiver  characteristic  values  are  unknown. 

2.  Two  options  in  the  type  of  interference  calculation 
desired.  One  is  strictly  a deterministic- type  calculation  between 
a given  transmitter  and  receiver.  The  other  is  a probabilistic- 
type  calculation  using  one  of  the  two  AVPAK  data  bases.  This 
calculation  yields  a graph  from  which  the  user  may  determine  the 
probability  of  interference  between  a transmitter-receiver  pair 
and  the  alterations  to  the  transmitted  power  which  must  be  made 

in  order  to  reduce  the  probability  of  interference. 

3.  The  development  of  two  data  bases  associated  with 
the  AVPAK  model,  viz.  AVBASE  and  AVFILE. 

a.  AVBASE  contains  nominal  characteristics  of 
selected  types  of  avionics  transmitters  and  receivers.  These 
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characteristics  may  be  used  in  running  the  AVPAK  model  for  a 
deterministic  type  run.  This  file  has  been  updated  to  include 
newly  manufactured  equipment  and  currently  contains  208  transmitters 
and  392  receivers.  The  equipment  information  in  AVBASE  is  also 
contained  in  the  ECAC  Nominal  Characteristics  File  (NCF) . A 
listing  of  the  ECAC  NCF  can  be  provided,  if  desired. 

b.  AVFILE  is  a functional  file  used  for  probabilistic 
processing.  In  AVFILE,  avionic  transmitters  and  receivers  are 
grouped  according  to  function.  The  group  is  represented  by  a 
single  symbolic  piece  of  equipment.  Its  characteristics  are  then 
determined  from  the  median  values  of  the  equipment  comprising  the 
group.  In  order  to  have  the  functional  file  as  statistically 
representative  as  possible  of  the  equipment  in  each  particular 
group,  work  has  been  done  to  redefine  certain  parameters  which 
are  used  to  represent  the  groups.  This  has  resulted  in  two  receiver 
groups  being  subclassified  according  to  their  spurious-response 
floor  levels. 


4.  A TSO/ARINC  characteristic  indicator  for  equipment 
located  in  AVBASE.  In  investigating  methods  to  make  the  AVFILE 
parameters  more  concise,  avionic  regulatory  standards  were  investi- 
gated with  the  hope  that  grouping  equipment  by  certain  standards 
would  yield  more  ’-epresentative  groupings.  This  did  not  prove 
advantageous.  As  a special  option,  the  user  may  inspect  a listing 
of  AVBASE  to  determine  what  equipment  has  associated  with  it  a 

TSO  (Technical  Standard  Order)  or  an  ARINC  (Aeronautical  Radio, 

Inc)  characteristic. 

5.  Calculation  of  power  density  on  or  near  the  airframe 
from  antennas  located  on  the  airframe  as  well  as  from  transmitters 
located  on  neighboring  aircraft.  Output  values  are  calculated 

in  both  dBm/m2  and  volts/m2  and,  if  appropriate,  the  model  will 
also  calculate  a cumulative  value  of  the  effects  of  more  than 
one  transmitter. 

6.  The  choice  of  three  coordinate  systems  for  user  input; 
(x,Y ,2) , (p ,d,z),  (B,w,F).  The  model  converts  the  latter  two  to 
cartesian  coordinates. 

7 . Validation  of  the  coupling  model . A total  of  996 
coupling  measurements  were  used  to  validate  the  model,  616  in 
the  UHF  portion  of  the  spectrum  and  380  in  the  SHF  portion.  For 
the  UHF  portion,  a comparison  of  measured  and  AVPAK- calculated 
data  yielded  a mean  error  of  -0.63  dB  with  a standard  deviation 
of  5.25  dB.  For  the  SHF  portion,  the  mean  error  was  '•■0.93  dB 
and  the  standard  deviation  was  5.7  dB.  For  the  combined  UHF  and 
SHF  data,  the  mean  error  was  -0.019  dB  with  a standard  deviation 
of  5.52  dB. 
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APPENDIX  A 

GEOMETRICAL  THEORY  OF  DIFFRACTION 


J.  B.  Keller  in  the  early  1950' s developed  the  "geometrical 
theory  of  diffraction"  ("GTD")  to  account  for  the  optical  phenomenon 
of  diffraction. 16  In  addition  to  the  usual  rays  of  geometrical 
optics,  GTD  introduces  new  rays  which  can  travel  along  curved  lines 
and  which  can  account  for  reflection  from  edges,  corners,  etc. 

These  new  rays  obey  several  laws  of  diffraction  which  are  analogous 
to  the  laws  of  reflection  and  refraction.  All  the  fundamental 
principles  of  ordinary  geometrical  optics  can  be  extended  to  geo- 
metrical diffraction.  The  only  difficulty  occurs  in  obtaining 
the  initial  value  of  the  field  at  the  point  of  diffraction.  For 
ordinary  rays,  the  field  of  a ray  emerging  from  a source  is  speci- 
fied at  the  source  but  for  a reflected  or  transmitted  ray,  the 
initial  value  is  obtained  by  multiplying  the  field  of  the  incident 
ray  by  a reflection  or  transmission  coefficient.  For  diffracted 
waves  the  initial  values  of  the  field  are  obtained  by  multiplying 
the  field  of  the  incident  rays  by  a diffraction  coefficient.  There 
are  different  diffraction  coefficients  for  edges,  vertices,  curved 
surfaces,  etc.  These  diffraction  coefficients  are  determined 
by  the  direction  of  incidence  and  diffraction,  the  wavelength,  and 
the  geometrical  and  physical  properties  of  the  medium  at  the  point 
of  diffraction  (see  Reference  4). 


In  an  ECAC  analysis  the  wing  was  modeled  as  a wedge.  Figure 
A-l  shows  the  dimensions,  angles  and  coordinate  systems  for  two 
antennas  on  opposite  sides  of  a wedge  of  angle  a.  In  Sach's 
development  of  the  diffraction  by  a single  wedge  (see  Reference 
4),  the  shadowing  attenuation  is  determined  by  the  following 
equation. 


A 


s 


where 


sin2  n/n 
n* 


X 


■(h  * hY  ♦ df 

i. 


(A-l) 


Ag  ■ shadowing  attenuation  of  a single  wedge  in  dB 


n 


variable  determined  from  the  wedge  angle,  a,  such 
that  n * --* 
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16Keller,  J.  B.,  The  Geometrical  Theory  of  Diffraction,  Symposium  on 
Microwave  Optics,  McGill  University,  Montreal,  Canada,  June,  1953. 
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A = transmitted  wavelength 

fcj  = minimum  distance  between  antenna  #1  and  apex  of 
wedge 

&2  = minimum  distance  between  antenna  #2  and  apex  of 

wedge 

d = distance  between  antenna  #1  and  antenna  #2  along 
z axis,  i.e.,  z2  - Zj 

Dq  = diffraction  coefficient 
and 


6 * angle  beyween  y axis  and  incident  ray 

0d  = angle  between  y axis  and  diffracted  ray  as  shown 
in  Figure  A-l . 

The  negative  sign  is  associated  with  soft  screen  boundary 
conditions,  i.e.,  field  component  tangent  with  respect  to  the 
diffracting  surface  and  the  positive  sign  is  associated  with  hard 
screen  boundary  conditions  i.e.,  the  field  is  normal  with  respect 
to  diffracting  surfaces. 

For  reasons  explained  in  Section  2,  DEVELOPMENT  OF  AN  AIRFOIL 
OBSTRUCTION  LOSS,  this  technique  did  not  meet  ECAC's  objective  and 
an  alternative  method  was  sought. 
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APPENDIX  B 

DESCRIPTION  OF  PARAMETERS  USED  TO  CREATE 
THE  FUNCTIONAL/PROBABILISTIC  FILE 


AVPAK  COUPLING  LOSS  ERROR  STATISTICS 

Six  sets  of  measurements  have  been  examined  to  determine  the 
error  to  be  expected  in  the  predicted  coupling  loss  determined 
with  the  automated  AVPAK  coupling- loss  routine.  Five  of  the  data 
sets  include  616  measurements  made  at  operating  frequencies  in 
the  UHF  portion  of  the  spectrum. 17> 18» 19,20,21 

Two  of  these  UHF  sets  constitute  a new  validation  of  the 
model;  the  remaining  UHF  and  SHF  sets  were  previously  used  in 
the  validation  of  the  AVPAK  2 model. 

The  sixth  set  of  data22  includes  380  measurements  made  at 
SHF  operating  frequencies  between  2.0  and  9.6  GHz.  These  measure- 
ments were  part  of  a classified  program  and  are  not  contained  in 
a published  report. 

UHF  Measurements  Results 

An  examination  of  this  data  indicates  that  the  mean  error 

(x)  is  -.63  dB  and  the  standard  deviation  of  the  error  is  5.25  dB. 
As  used  herein,  a negative  error  indicates  that  the  measured 
coupling  loss  was  lower  than  the  AVPAK  predicted  coupling  loss 
and  a positive  error  indicates  that  the  measured  loss  was  higher 
then  the  predicted  value. 


17Electronic  Communications,  Inc.,  Electromagnetic  Compatibility 
Report  for  KC-135B  Aircraft,  January  1965. 

10The  Boeing  Corporation,  Category  II  Flight  Test  Report  for 
KC-135B  (PACCS)  Electronic  System,  Test  No.  T6-3181,  1965. 

19Mart  in,  H.,  Measured  Adjacent  Signal  Interference  of  Collocated 
AN/ ARC-51  Transceivers  (U) , ESD-TR-67-003,  January  1968, 
CONFIDENTIAL. 

20The  Boeing  Corporation,  EC-135C  AFSAT  EMC  Baseline  Measurements 
and  Analysis,  Test  No.  T3-1702,  July  1974. 

21E-System  Inc.,  Garland  Division,  Model  No.  E-4A,  Report  No. 
G8494. 12.26,  1973. 

22Zimballatti,  A.,  Grumman  Aircraft  Corporation,  Personal  Contact, 
June  1972. 
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The  error  distribution  appears  to  be  normal  with  a median 
error  of  0 dB  and  a standard  error  of  6.0  dB,  based  on  a x2 
(chi-square)  test  with  6 degrees  of  freedom  at  a .05  signifi- 
cance level.  The  range  of  errors  observed  is  shown  in  Figure 
B-l . 

SHF  Measurement  Results 

The  SHF  data  fell  into  two  categories.  The  first  category 
included  292  measurements  with  antennas  situated  such  that  the 
measurements  provided  comparisons  of  observations  versus  pre- 
dictions from  the  AVPAK  knife-edge  coupling  subroutine.  The 
remaining  88  measurements  were  made  between  antenna  pairs  with 
both  antennas  mounted  aft  of  the  nose  bulkhead.  The  second 
category  of  data  enabled  an  examination  of  the  raised  antenna 
versus  surface-mounted  antenna  subroutine  and  the  routine  which 
calculates  coupling  loss  along  a curved  surface. 

The  range  of  errors  observed  for  the  first  set  of  SHF  data 
is  shown  in  Figure  B-2. 

The  results  of  the  first  category  of  SHF  measurements  indicate 

that  the  mean  error  (x)  is  +1.05  dB;  the  standard  deviation  of  the 
error  (a)  is  6.0  dB.  The  error  distribution  was  tested  using  a 
X2  (chi-square)  test  at  a .05  significance  level  using  6 degrees 
of  freedom.  The  results  indicate  that  these  errors  are  normally 
distributed  with  a median  error  of  0 dB  and  a standard  deviation 
of  6.0  dB. 

The  mean  error  (x)  of  the  second  set  of  SHF  measurements  was 
found  to  be  +.07  dB  and  the  standard  deviation,  4.9  dB. 

When  both  SHF  categories  were  considered  as  one  group,  it 

was  found  that  the  mean  error,  (x)  was  +0.93  dB  and  the  standard 
error  (o),  5.7  dB.  The  distribution  of  the  combined  SHF  errors 
appears  to  be  normal  with  a median  error  of  0 dB  and  a standard 
error  of  6.0  dB,  based  on  a x2  (chi-square)  test  at  .05  signifi- 
cance level  and  6 degrees  of  freedom.  The  range  of  errors  for 
the  combined  380  SHF  measurements  may  be  observed  in  Figure  B-3. 

SUWlARy  OF  COMBINED  RESULTS 

A total  of  996  measurement  points  were  evaluated  to  determine 
the  error  to  be  expected  in  the  AVPAK  coupling  loss  predictions. 

The  range  of  observed  errors  of  all  the  measurements  is  presented 
in  figure  B-4.  The  distribution  of  the  expected  errors  appears  to 
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Figure  B-l.  Observed  errors  in  UHF  coupling  loss  prediction. 


Figure  B-3.  Observed  errors  of  predicted-versus-measured  coupling  loss  in  the 
SHF  band  - all  paths. 


Figure  B-4.  Observed  errors  in  AVPAK  coupling  loss  prediction, 
coabined  results. 
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APPENDIX  C 

GENERAL  FLOW  CHART  FOR  AVPAK  3 
(Page  1 of  3) 
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APPENDIX  D 

LISTING  OF  THE  CONTENTS  OF  AVBASE  AND  AVFILE 


The  information  in  the  AVBASE  and  AVFILE  files  is  intended  to 
be  used  as  input  data  to  the  AVPAK  3 model.  It  may  not  be  directly 
comparable  to  published  specifications.  Many  parameters  were 
synthesized  using  equipment  schematics,  for  reasons  discussed  on 
pages  10  and  11.  Abbreviations  used  in  the  AVBASE  file,  listed 
first,  are  defined  in  TABLE  D-l. 

In  the  computer  printout  that  follows  TABLE  D-l,  entries 
of  ".0,"  ".00"  and  ".000"  should  be  read  as  blanks. 
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BW1 : 

BN2: 

EBW: 

FLD: 

FP: 

HI  FLT: 

HI  FRQ: 

HI  SPR  FRQ: 

IF: 

IF  BN: 

IF  SF: 

IM  REJ: 

INT: 


TABLE  D-l 

AVBASE  ABBREVIATIONS 
(Page  1 of  4) 


The  emission  envelope  bandwidth  of  the  first  break- 
points. This  is  the  nominal  3 dB  bandwidth  in  MHz 
(pp.  35-36). 

The  emission  envelope  bandwidth  of  the  second 
breakpoints.  This  is  the  bandwidth  at  which  the 
envelope  shows  a second  fall-off  characteristic 
and  is  in  MHz. 

Emission  designator  bandwidth.  It  is  the  bandwidth 
containing  99%  of  the  mean  radiated  power.  Its 
units  are  MHz. 

Field  format. 

Floating  point  format. 

Upper  filter  frequency  in  MHz.  The  lower  and  upper 
filter  limits  are  used  to  truncate  the  spectral 
energy  of  the  transmitted  emissions  outside  the 
two  frequencies  given.  The  primary  application 
for  these  data  is  to  describe  waveguide  cut-off 
phenomena. 

Upper  operating  frequency  of  the  equipment.  Its 
units  are  MHz. 

Upper  spurious  response  limit.  Similar  to  LO  SPR 
FRQ  but  related  to  upper  side  of  the  RF  selectivity 
curve.  Its  units  are  in  Miz. 

Intermediate  frequency.  It  is  the  first  intermediate 
frequency  if  there  is  more  than  one.  Its  units  are 
MHz. 

Intermediate  frequency  bandwidth.  This  is  the  final 
IF  3 dB  bandwidth.  Its  units  are  MHz. 

Intermediate  frequency  (final  IF)  selectivity  slope. 
Its  units  are  dB/decade.  (See  IF  Skirt  Slope,  p.  29.) 

Image  rejection.  Its  units  are  dB  (p.  33). 

Integer  format. 
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KEY  NO: 

LO  FLT: 

LO  FRQ: 

LO  POS: 

LO  SPR  FRQ: 

MOD  TYP: 


TABLE  D-l 
(Page  2 of  4) 


The  key  number  of  a unique  equipment  I.D.  assigned 
in  alpha-numeric  order  according  to  manufacturer 
and  equipment  function  (pp.  54-55). 

Lower  filter  frequency  in  MHz  (see  HI  FLT,  on 
previous  page) . 

Lower  operating  frequency  of  the  equipment.  Its 
units  are  MHz. 


Local  oscillator  position.  A = above,  B = below, 
relative  to  the  carrier  frequency.  If  an  input  is 
not  given  it  is  assigned  a "C"  and  treated  as  being 
both  above  and  below  the  carrier. 

Lower  spurious  response  limit.  It  is  a discrete 
frequency  defined  by  the  intersection  of  the  spurious 
floor  with  the  lower  frequency  side  of  the  RF  selec- 
tivity curve.  Its  units  are  in  MHz  (P.  35). 


FCC  emission  modulation  type.  Emissions  are  classified 
and  symbolized  according  to  the  following  character- 
istics: 


Types  of  modulation  of  main  carrier: 


Symbol 


a)  Amplitude  A 

b)  Frequency  (or  Phase)  F 

c)  Pulse  P 

Types  of  transmissions: 


a)  Absense  of  any  modulation  intended  to  0 

carry  information. 

b)  Telegraphy  without  the  use  of  a modu-  1 

lating  audio  frequency. 

c)  Telegraphy  by  the  on-off  keying  of  a 2 

modulating  audio  frequency  or  audio  fre- 
quencies, or  by  the  on-off  keying  of  the 
modulated  emission.  (Special  case:  an 

unkeyed  modulated  emission. ) 

d)  Telephony  (including  sound  broadcasting) . 3 

e)  Facsimile  (with  modulation  of  main  carrier  4 

either  directly  or  by  a frequency  modu- 
lated sub-carrier). 
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(Page  3 of  4) 


Symbol 

f)  Television  (vision  only)  5 

g)  Four- frequency  diplex  telegraphy  6 

h)  Multichannel  voice- frequency  telegraphy  7 

i)  Cases  not  covered  by  the  above  9 


Supplementary  characteristics: 


a) 

Double  sideband 

(none) 

b) 

Single  sideband: 

reduced  carrier 

A 

full  carrier 

H 

suppressed  carrier 

J 

c) 

Two  independent  sidebands 

B 

d) 

Vestigial  sideband 

C 

e) 

Pulse: 

amplitude  modulated 

D 

width  (or  duration)  modulated 

E 

phase  (or  position)  modulated 

F 

code  modulated 

G 

NOMENCLATURE:  1)  First  part  of  column:  Manufactures 

unique  name  of  equipment,  further 
identified  by  "KEY  NO."  to  manufacturer 
(see  p.  55) . 

2)  Second  part  of  column:  Type  of  equip- 

ment (i.e.,  C = VHF  Comm.;  VL  = VOR, 
Localizers;  G = Glideslope;  ATC  = Air 
Traffic  Control  Transponders;  WR  = 
Weather  Radar;  DME  = Distance  Measuring 
Equipment;  ALT  = Altimeter). 

PCI:  Pulse  compression  indicator.  This  requires 

a "C"  for  chirp  modulation  only,  otherwise 
this  field  is  blank. 

PRF:  Pulse  repetition  frequency.  Its  units  are 

Hz. 

PW:  Pulse  width.  It  is  the  value  between  half 

amplitude  points.  For  chirped  modulation 
it  is  the  stretched  width.  Its  units  are 
us.  (See  p.  38  for  explanation  of  stretched 
width. ) 
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(Page  4 of  4) 


PWR: 


R/F: 


RF  SF: 
SENS: 

SF1 : 

SF2: 

S/I: 

SP  REJ: 
TSO  IND: 


The  transmitter  power.  It  is  the  average  power 
output  for  communication  transmitters  and  the 
peak  power  output  for  pulsed  transmitters.  Its 
units  are  dBm. 

The  average  pulse  rise  and  fall  time  (6).  It  is 
determined  from: 


6 


2 

l/«r  ♦ l/«f 


where  6 = the  time  for  the 

r 


pulse  to  rise  from  its  10%  amplitude  to  its  90% 
value  and  6^  = the  time  to  fall  from  90%  to  10%. 

Its  units  are  us. 


RF  selectivity  slope.  Its  units  are  dB/decade 
(p.  29). 

Receiver  sensitivity.  It  must  be  preceded  by  a 
negative  sign.  Its  units  are  aBm. 

The  first  slope  fall-off  of  the  emission  envelope. 
It  is  the  slope  adjacent  to  the  3 dB  bandwidth. 

Its  units  are  dB/decade  (see  p.  37). 

The  second  slope  fall-off  of  the  emission  envelope. 
It  is  the  slope  at  frequencies  greatly  separated 
from  the  tuned  frequency.  Its  units  are  dB/decade 
(see  M^,  p.  37) . 

Degradation  threshold.  It  is  the  minimum  signal- 
to- interference  ratio  required  for  non-interference 
Its  units  are  dB  (p.  43). 

Spurious  response  rejection.  Its  units  are  dB 
(P-  34). 


Technical  Standard  Order  indicator.  T = equipment 
is  TSO'd.,  A = equipment  is  TSO'd  and  also  falls 
under  an  ARINC  characteristic  (see  pp.  10  and  11). 
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APPENDIX  E 

ACTUAL  PROGRAM  RUN  DECKS  AND  OUTPUT 


INTRODUCTION 

There  are  twenty- four  possible  types  of  AVPAK  program  execu- 
tions, as  shown  in  TABLE  E-l.  Each  type  of  nan  is  different  with 
respect  to  the  type  of  calculation  desired  (power  density  or  INR), 
the  type  of  analysis  desired  (cosite  or  intersite),  the  type  of 
INR  answer  desired  (deterministic,  functional  or  probabilistic), 
or  the  coordinate  system  used  for  the  inputs,  (aircraft  industry, 
cylindrical,  or  rectangular).  Similarly,  there  are  twelve  impos- 
sible run  type  combinations,  as  shown  in  TABLE  E-2. 

SAMPLE  RUN  DECKS  AND  OUTPUT 

Seven  different  AVPAK  sample  runs  are  discussed  below.  Each 
run  corresponds  to  one  or  more  of  the  twenty-four  possible  types. 
Together  they  exhibit  extensive  execution  of  the  capabilities  of 
the  AVPAK  model,  and  they  provide  the  user  with  an  illustration 
of  actual  sets  of  data  cards.  The  output  corresponding  to  each 
run  is  also  shown. 

SAMPLE  RUN  #1  (TYPE  #17  ANALYSIS) a 

This  sample  analysis  involves  one  transmitter  with  a direc- 
tional antenna  on  an  aircraft  and  one  receiver  with  an  omnidirec- 
tional antenna  on  the  ground.  The  input  data  cards  are  shown  in 
Figure  E-l,  and  describe  th^  situation  that  follows. 

General  Parameter  Data  * 

This  run  is  to  perform  an  intersite,  deterministic,  inter- 
ference-to-noise  ratio  analysis.  Inputs  are  in  the  rectangular, 
(X,Y,Z),  coordinate  system.  The  radius  of  the  fuselage,  although 
not  used  in  an  intersite  analysis,  is  given  as  100  inches. 

Intersite  Data 

Site  1 heading  « 50° 

Site  1 altitude  * 11,560  feet 

Site  2 heading  = 08 

Site  2 altitude  * 1.0  feet 


aRefer  to  TABLE  E-l  for  description  of  types  of  runs. 
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TABLE  E-l 


ALLOWABLE  RUN  TYPES 


Type 

Calculation 

Analysis 

Answer 

Coordinate  System 

1 

Power  Density 

Cosite 

N/Aa 

(B.W.F) 

2 

M If 

II 

II 

(X.Y.Z) 

3 

M II 

♦ 1 

II 

(p.e.z) 

4 

II  II 

Intersite 

II 

(B.W.F) 

5 

If  II 

II 

If 

(X.Y.Z) 

6 

If  If 

II 

II 

(p.e.z) 

7 

INR 

Cosite 

Deterministic 

(B.W.F) 

8 

II 

II 

If 

(X.Y.Z) 

9 

II 

II 

II 

(p.e.z) 

10 

ft 

II 

Functional 

(B.W.F) 

11 

II 

If 

II 

(X.Y.Z) 

12 

II 

ft 

It 

(p.e.z) 

13 

II 

II 

Probabilistic 

(B.W.F) 

14 

II 

II 

II 

(X.Y.Z) 

IS 

II 

II 

II 

(p.e.z) 

16 

II 

Intersite 

Deterministic 

(B.W.F) 

17 

II 

II 

II 

(X.Y.Z) 

18 

II 

II 

II 

(p.e.z) 

19 

II 

If 

Functional 

(B.W.F) 

20 

II 

II 

II 

(X.Y.Z) 

21 

II 

II 

II 

(p.e.z) 

22 

ft 

II 

Probabilistic 

(B.W.F) 

23 

11 

II 

II 

(X.Y.Z) 

24 

II 

II 

I# 

(p.e.z) 

“The  answer-type  terminology  applies  only  to  INR  runs.  With  Power 
Density  runs  the  terminology  is  not  applicable. 


TABLE  E-2 


IMPOSSIBLE  RUN  TYPES 


Calculation 

Analysis 

answer 

Coordinate  System 

Power  Density 

Cosite 

Functional 

(B.W.F) 

II  II 

II 

II 

(X.Y.Z) 

I ••  II 

If 

(p.e.z) 

II  II 

II 

Probabilistic 

(B.W.F) 

If  II 

II 

II 

(X.Y.Z) 

•I  II 

II 

II 

(P.e.z) 

• 1 II 

Inters ite 

Functional 

(B.W.F) 

II  If 

It 

II 

(X.Y.Z) 

II  II 

II 

II 

(P.e.z) 

II  II 

II 

Probabilistic 

(B.W.F) 

• 1 II 

If 

•I 

(X.Y.Z) 

II 

II 

(p.e.z) 
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Intersite  Data  (Cont'd.) 

Ground  distance  from  site  1 to  site 
True  bearing  from  site  1 to  site  2 
Platform  indicators:  site  1 

site  2 


2 = 10.0  statute  miles 
= 050° 

= moving  (airplane); 
= fixed  (ground) 
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The  antenna  for  site  number  2 is  a vertically  polarized  di- 
pole and  is  considered  to  be  omnidrectional . 

Transmitter  Data 


The  single  set  of  two  transmitter  data  cards  describe  the 
characteristics  of  the  equipment  named  WEATHER  RADAR.  The  second 
data  card  of  this  set  is  required  to  enter  the  transmitter  filter 
frequency  limits. 

Receiver  Data 


The  single  receiver  data  card  describes  the  characteristics 
of  the  equipment  named  DOPPLER  RADAR. 

The  data  of  Figure  E-l  have  been  applied  to  the  diagram  in 
Figure  E-2,  to  demonstrate  a portion  of  the  logic  path  taken  for 
the  determination  of  the  off-axis  angles,  required  to  calculate 
the  transmitter  antenna  gain.  Figure  E-3  is  an  illustration  of 
the  intersite  geometry  involved  in  the  processing,  and  Figure  E-4 
contains  the  actual  printout  results  of  the  executed  program.  As 
can  be  seen  from  the  output,  the  following  calculated  or  input 
values  were  used  in  the  EMC  analysis: 

PT  = transmitter  power  (input)  = 73.0  dBm 

GT  = transmitter  antenna  gain  (calculated)  = 15.6  dBi 

GR  = receiver  antenna  gain  (assigned)  = 2.0  dBi 

LP  = site-to-site  path  loss  (calculated)  = 130.3  dB 

LF  = transmitter-receiver  off-frequency  rejection 
(calculated)  = 80.0  dB 

S/I  = signal-to-interference  ratio  threshold  (input) 

= 10.0  dB 

RS  = receiver  sensitivity  level  (input)  = -120.0  dBm 

PI  = interfering  power  level  at  the  receiver  input 
terminals  (calculated) 

= -140.8  dB. 

It  can  be  seen  that  two  receiver  characteristic  values  were 
reset  or  recovered  by  the  program:  an  RF  slope  fall-off  of  200 

dB/decade  and  an  image  rejection  of  60  dB. 
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( MAMMOMAM  • CMC  ANALYSIS) 


Figure  E-2.  Determination  of  the  difference  angles  6 and  0, 

between  the  transmitter  antenna  main  beam  and  the 
transmitter-site  to  receiver-site  path. 
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SAMPLE  RUN  #2  (TYPE  #1  RUN) 


This  sample  analysis  involves  six  transmitters  and  a single 
power  density  point,  all  on  the  same  aircraft.  The  input  data 
cards  are  shown  in  Figure  E-5,  and  serve  to  describe  an  analysis 
that  follows. 

General  Parameter  Data 

This  run  is  to  calculate  the  power  density  at  a point  on  an 
aircraft,  resulting  from  six  collocated  transmitters.  The  inputs 
are  with  respect  to  the  aircraft  industry's  butt-line,  water  line, 
fuselage  station,  (B,9f,F),  coordinate  system.  The  maximum  fuse- 
lage radius  is  140.32  inches.  No  bulkhead  obstruction  on  the 
fuselage  is  assumed  to  exist. 


{weapon!  l 

lEMD  1 1 
I i Hirx 

I 2 Hf* 

I 3 ll  (TX 


7100  £140.32  SAMPLE  RUM  «2  (TYPE  1) 

1 195.  133.16  382.2  ZORRU  PORT  1 4.8 

X.  AMT.  1 4 6.0  0.  208. < 

X.  AMT.  2 4 6.0  0 . 208.1 

X.  AMT.  3 4 6.9  0.  203.< 

X.  AMT.  4 2 6.0  0.  280.1 

X.  AMT.  5 4 6.0  0.  120.: 

X.  AMT.  6 4 -5.  9.  120/ 


6 1|  jTX.  ANT.  6 

7 l|  MEAPOM  PT. 
1APN-25  A 
2AFft-25  B 

s 

3Af£t-34  C 
6ARH-21 
7HERPOM  PT. 


300.0 

400.0 

255.0 


6.0 

0.  208.6 

32.2 

6.0 

0.  208.6 

80.5  | 

6.9 

0.  203.6 

333. 0 | 

6.0 

0.  280.6 

545.0  | 

6.0 

0.  120.7 

557.5  | 

-5. 

0.  120.7 

163.0  i 

4.8  4.8 

6.0  i 

1 

1090.0 

60. 

R0  1.0  | 

14. 

1090.0 

60. 

P0  1.0  | 

14. 

300.0 

60. 

A1  i 

0. 

400.0 

60. 

A1  I 

255.0 

60. 

A1  | 

1025.0 

60. 

P0  3.5  i 

.15 

FORTRAN  STATEMENT 
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Pod  Obstruction  Data 


One  pod  is  to  be  considered  in  the  analysis,  ZORRO  PORT  1. 

The  nose-centroid  of  this  pod  is  located  at  coordinates  (195 
inches,  133.16  inches,  382.2  inches)  in  the  (B,W,F)  coordinate 
system  of  the  fuselage.  The  pod  is  assumed  to  be  a cylinder  of 
radius  4.8  inches.  / 

i 

Antenna  and  Power  Density  Point  Pat* 

c 

Each  transmitter  is  coupled  to  its  own  antenna.  The  antenna 
types  are  all  undefined.  One  is  mounted  on  the  vertical  stabilizer, 
and  the  others  are  mounted  on  the . fuselage. . One  antenna  has  an 
input  gain  of  -5  dBi,  the  others  have  an  input  gain  of  6 dBi. 

The  location  of  each  antenna  is  given  in  the  (B,W,F)  coordinate 
sytem  of  the  fuselage.  The  location  of  tha  power  density  point 
is  specified  with  respect  to  the  pod  in  the  pod's  (B,Hf,F)  coordinate 
system. 

Transmitter  Data 

The  nomenclatures  and  cha^acfcariftics  of  the  six  transmitters 
are  given  on  six  cards.  Each  transmitter  is  associated  with  a 
different  antenna  by  use  of  a unique  antenna  identification  number. 
Each  transmitter  frequency,  power,  and  modulation  type  is  specified. 
The  peak  transmitter  power  is  60  dBm  for  all  equipments.  The 
pulsewidth,  pulse  rise/fall  times,  and  pulse  repetition  frequency 
for  each  of  the  pulsed  (modulation  type  PO)  transmitters  are  given. 

The  power  density  point  data  card  specifies  the  nomenclature 
and  point  identification  number  Note  that  the  power  density 
point  is  treated  the  same  as  an  antenna  to  simplify  the  corre- 
spondences of  point  nomenclature  and  point  location  coordinates. 
Obviously  the  antenna  and  equipment  characteristic  data  are  not 
applicable  for  the  power  density  point,  so  that  most  data  fields 
on  both  cards  are  blank. 

Figures  E-6  and  E-7  illustrate  the  geometry  of  the  power 
density  analysis  to  be  performed.  Figure  E-8  shows  the  results 
of  the  power  density  analysis.  The  expected  peak  and  average 
power  densities  at  the  pod  due  to  each  transmitter  individually, 
and  then  the  cumulative  power  due  to  all  transmitters  are  shown. 

Note  that  the  peak  power  density  computation  is  not  appropriate 
for  non-pulsed  equipments.  The  field  strengths  in  volts/meter 
corresponding  to  each  power  density  level  in  dBm/meter2  are  also 
listed.  Note  that  messages  produced  by  the  model  indicate  that 
buttline  coordinates  were  not  given  for  the  six  antenna  locations. 
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Figure  E-7.  Location  of  the  power  density  point  on  the  Pod  "ZORRO  PORT  1”  in  Sample 


Figure  E-8.  Printout  of  Sample  Run  #2. 
(Page  1 of  2) 
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These  particular  messages  are  diagnostics  for  the  user's  benefit 
but  should  not  be  of  any  concern  in  this  analysis,  as  the  butt- 
line dimensions  of  all  antennas  were  intended  to  be  zero. 


SAMPLE  RUN  #3  (TYPE  #9  RUN 


The  data  cards  for  this  sample  run,  shown  in  Figure  E-9 
describe  an  analysis  to  be  performed  that  follows. 


ISO.  45.6  SAMPLE  HUM  *3  < TYPE  9t 
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General  Parameters 


This  run  is  to  perform  a cosite,  deterministic,  INR  analysis 
of  five  transmitters  and  two  receivers.  Input  units  are  to  be  in 
the  cylindrical  (p,  0,  Z)  coordinate  system.  The  maximum  fuse- 
lage radius  is  152.4  inches.  A bulkhead  obstruction  with  a height 
of  45.6  inches  exists  150.0  inches  aft  of  the  fuselage  nose.  There 
are  no  wing  or  pod  obstructions  to  be  considered  in  the  analysis. 

Antenna  Data 


Each  equipment  uses  a different  antenna.  All  antennas  are 
vertically  polarised,  fuselage-mounted  dipoles.  The  locations  of 
the  seven  antennas  are  specified  in  the  p,  0,  Z coordinate  system. 
The  seven  antennas  are  numbered:  1,  2,  3,  4,  5,  7,  and  8. 

Transmitter  Data 


The  five  transmitters  are  tuned  to  225,  1025,  or  1090  MHz. 
All  transmitters  are  assumed  to  have  the  same  first  and  second 
emission  bandwidth,  power,  and  number  of  harmonics  to  be  con- 
sidered in  the  analysis.  All  transmitters  except  the  IFF  are 
assumed  to  have  the  same  falloff  characteristics  and  modulation 
types. 

Receiver  Data 


One  receiver  is  tuned  to  118  MHz  and  the  other  to  120  MHz. 

The  local  oscillator  tracks  above  the  tuned  frequency  for  both 
receivers.  All  other  receiver  characteristics  are  to  be  retrieved 
from  the  AVBASE  file  from  those  records  with  identification  num- 
bers 174  and  1752. 

The  results  of  the  INR  analysis  are  shown  in  Figure  E-10. 

Note  that  no  transmitter-receiver  combinations  are  potential 
interference  problems,  as  all  ten  combinations  are  assumed  to  be 
interference- free.  All  values  for  the  terms  in  the  interference 
expression  are  calculated  as  being  equal  except  for  the  path 
losses  and  off-frequency  rejections.  The  first  harmonic  (funda- 
mental) was  determined  to  be  the  most  critical  transmitting  fre- 
quency to  be  considered  in  calculating  the  off-frequency  rejections. 

The  locations  of  the  antennas  on  the  fuselage  are  shown  in 
the  diagram  of  Figure  E-U,  along  with  the  calculated  path  losses. 

As  in  Sample  Run  #2,  several  diagnostic  messages  pertaining  to 
antenna  location  data  are  printed.  Also  listed  are  the  transmitter 
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Figure  E-10.  Printout  of  Sample  Run  #3. 
(Page  1 of  3) 
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spectral  slope  fall-off  values  recovered  in  the  program.  The 
receiver  characteristics  as  retrieved  from  AVBASE  are  listed  under 
the  "Receivers  with  Good  Data." 

SAMPLE  RUN  *4  (TYPE  »15  RUN) 

The  data  cards  for  this  sample  run,  shown  in  Figure  E-12, 
describe  an  analysis  to  be  performed  that  follows. 
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Figure  E-12.  Data  card  deck  for  Sample  Run  #4. 


General  Parameters 

This  sample  is  for  a cosite,  probabilistic,  INR  analysis. 

In  this  run  one  transmitter  and  one  receiver  are  involved.  The 
input  units  are  in  the  cylindrical  coordinate  system.  A fuselage 
of  56  inches  radius  with  a 16.56-inch  bulkhead  located  at  36 
inches  aft  of  the  fuselage  nose  is  to  be  considered.  No  other 
obstructions  are  to  be  considered  in  the  analysis. 

Antenna  Data 

The  transmitter  antenna  ia  a horizontally  polarized  parabolic 
antenna  18  inches  in  diameter  located  near  the  nose  of  the  air- 
craft, forward  of  the  bulkhead  obstruction.  The  receiver  antenna 
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is  a horizontally  polarized  horn  antenna,  with  an  aperture  of 
4.3  inches  by  3.1  inches,  and  it  is  mounted  on  the  fuselage  aft 
of  the  bulkhead  obstruction. 

Transmitter  Data 


The  transmitter  uses  antenna  number  1.  The  equipment  char- 
acteristics are  to  be  retrieved  from  the  AVFILE  record  whose 
identification  number  is  5. 

Receiver  Data 

The  receiver  uses  antenna  number  2.  The  equipment  charac- 
teristics are  to  be  retrieved  from  the  AVFILE  record  whose  identi- 
fication number  is  57. 

The  locations  of  the  two  antennas  on  the  fuselage  with  re- 
spect to  the  bulkhead  are  shown  in  Figure  E-13.  The  output  from 
Sample  Run  #4  is  shown  in  Figure  E-14.  The  transmitter  and  re- 
ceiver characteristics  retrieved  from  AVFILE  are  listed. 

Because  the  receiver  input  power  level  of  -213.0  dBm  does 
not  exceed  the  receiver  sensitivity  of  -72.9  dBm,  the  predicted 
probability  of  interference  vs.  dB  error  is  not  plotted.  The 
table  of  interfering  equipments  shows  no  potential  interference 
between  the  two  equipments. 

SAMPLE  RUN  »5  (TYPE  »14  RUN) 

The  data  cards  for  this  sample  run,  shown  in  Figure  E-15, 
describe  the  following  analysis. 

General  Parameters 

Similar  to  Sample  Run  #4,  this  run  is  a cosite,  probabilistic, 
INR  analysis.  In  this  analysis,  however,  the  inputs  are  in  the 
rectangular,  (X,Y,Z),  coordinate  system,  and  no  bulkhead  obstruc- 
tion is  to  be  considered.  The  fuselage  radius  is  240  inches.  No 
wings  or  pods  are  to  be  considered. 

Antenna  Data 

The  transmitter  antenna  is  a vertically  polarized,  parabolic 
antenna  with  a 12 -inch  aperture,  mounted  on  the  fuselage.  The 
receiving  antenna  is  an  unspecified  type,  mounted  on  the  fuselage 
with  a gain  value  of  3 dBi.  The  type  is  unspecified  since  a gain 
value  is  entered. 
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Figure  E-15.  Data  cards  for  Sample  Run  #5. 


The  transmitter  uses  antenna  number  1 and  has  the  charac- 
teristics of  the  AVFILE  record  with  identification  number  2.  The 
receiver  uses  antenna  number  2 and  has  the  characteristics  of  the 
AVFILE  record  with  identification  number  51. 

The  output  for  Sample  Run  #5  is  shown  in  Figure  E-16.  Since 
the  effective  median  interference  power  of  -104.2  dBm  and  the 
upper  decile  cull  correction  factor  both  exceed  the  receiver 
sensitivity,  the  predicted  probability  of  interference  vs.  dB 
error  is  plotted.  The  actual  coordinates  of  21  points  on  the 
graph  are  also  printed.  The  table  of  interfering  equipments 
notes  that  these  two  equipments  are  a potential  source  of  inter- 
ference. 

SAMPLE  RUN  #6  (TYPE  #3  RUN) 

The  data  card  deck  for  this  sample  run  is  shown  in  Figure 
E-17.  This  run  performs  a cosite,  power  density  analysis  among 
five  transmitters  and  four  power  density  points.  Inputs  are  to 
be  in  the  cylindrical,  (p,  8,  Z),  coordinate  system.  The  maximum 
fuselage  radius  is  100  inches. 
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Figure  E-16.  Printout  of  Sample 
(page  1 of  3) 


FAA-RD-76-50 


Appendix  E 


Figure  E-16.  (Page  2 of  3) 
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Figure  E-17.  Data  card  deck  for  Sample  Run  #6. 


A bulkhead  of  height  40  inches,  located  50  inches  aft  of  the 
fuselage  nose,  is  to  be  considered.  An  airfoil  is  to  be  consid- 
ered as  a potential  obstruction  to  the  transmission  path.  The 
front  wing  edges  are  located  280  inches  aft  of  the  fuselage  nose. 
The  aft  wing  edges  are  located  400  inches  aft  of  the  nose.  The 
port  wing- fuselage  intersection  points  are  located  at  100°  with 
respect  to  the  vertical  stabilizer.  The  five  transmitter  antennas 
are  all  different  types  and  mounted  at  various  locations  about 
the  aircraft.  They  also  represent  a mixture  of  polarizations  and 
sizes.  The  four  power  density  points  are  also  located  at  dif- 
ferent points  about  the  aircraft. 
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Figure  E-18.  Location  of  transmitter  antennas  and  power  density  points  for  Sample  Run  #6 
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Figure  E-20.  Data  card  deck  for  Sample  Run  #7. 
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Figure  E-21.  Printout  of  Sanple  Run  #7. 
(Page  1 of  3) 
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